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Introduction 


1.  Title  of  the  Project 

Space  Weather  Effects:  Decay  Phases  in  Gradual  and  Impulsive  Solar 
Energetic  Particle  Events. 

2.  Contracting  Institute 

D.V.Skobeltsyn  Institute  of  Nuclear  Physics  of  Moscow  State  University 

Address:  1 19992,  Russia,  Moscow,  Leninskiye  gory 

Director,  Prof 

Panasyuk  Michail  Igorevich 

Telephone:  7-095-9391818 

Fax:  7-095-9390896 

E-mail;  mip@srdmail.sinp.msu.ru 

3.  Participating  Institutes 

No 

4.  Project  Manage 

Daibog  Elena  Isaevna,  Senior  Scientific  Researcher,  Ph.D 
Telephone;  7-095-9392548 
Fax:  7-095-9390896 

E-mail:  daibog@srd.sinp.msu.ru 

5.  Commencement  Date,  Duration 

June  1,  2001,  duration  -  24  months 


The  purpose  of  the  project  is  the  investigation  of  solar  energetic  particle  (SEP)  events 
with  the  practical  application  to  the  problems  of  space  weather  and  monitoring  of 
radiation  conditions  in  interplanetary  space.  Scientific  goals  of  the  project  are  related  to 
the  basic  problem  of  solar  cosmic  ray  physics  which  includes  determination  of  flare  vs. 
CME-  driven  shock  contributions  to  particle  fluxes  in  SEP  events  in  interplanetary  space 
(accelerated  in  flares  and  by  CME-driven  shocks),  forecasting  of  radiation-dangerous 
events  and  determination  of  properties  of  interplanetary  medium  (populated  by  the 
particles  observed)  by  the  parameters  of  particle  fluxes. 


Expected  results 

At  the  project  elaboration: 

•  the  .statistically  provided  data  base  of  SEP  events 
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•  definition  of  empirical  dependencies  of  SEP  event  decay  phase  characteristics 
(such  as  particle  flux  decay  rate,  energetic  particle  spectrum,  anisotropy  of 
particle  fluxes  etc.)  upon  a  distance  from  the  Sun  and  heliolongitude  under 
various  conditions  in  solar  corona  and  in  interplanetary  space; 

•  influence  of  flare  and  shock  accelerated  particles  on  the  observed  fluxes  at  decay 
phases  of  events; 

•  determination  of  conditions  of  spatial  and  temporal  invariance  of  particle  spectra 
during  SEP  event  decay  phase; 

•  detailed  examination  of  SEP  event  characteristic  decay  time  intensity  variations 
over  the  solar  activity  cycle 

From  the  practical  point  of  view  obtained  results  can  be  applied  to  space  weather  models 
if  someone  considers  space  weather  as  a  set  of  heliospheric  phenomena  including 
interplanetary  shocks,  recurrent  high  velocity  streams  in  solar  wind,  plasma  clouds  etc. 
In  particular,  SEPs  data  are  suitable  for  estimations  of  characteristic  spatial  scale  of 
interplanetary  shocks  and  their  parameter  (velocity  etc.)  variations  with  radial  and 
azimuthal  distance. 

Technical  approach 

The  Project  was  elaborated  using  the  whole  totality  of  the  data  sets  on  solar  wind, 
interplanetary  magnetic  field  and  solar  energetic  particles  that  was  obtained  during  long 
period  of  observations  in  space,  as  well  as  the  data  on  coronal  mass  ejections,  flares  and 
other  phenomena  of  solar  activity.  The  Project  foundation  is  the  data  of  m/s 
measurements  in  which  the  same  SEP  events  together  with  solar  wind  conditions  were 
observed  at  different  space  regions  simultaneously.  The  database  includes  the  results  of 
measurements  onboard  the  Helios- 1,2,  IMP-5,7,8,  GRANAT  and  ISEE-3,  as  well  as 
data  of  recent  observations  onboard  SOHO  and  WIND  s/c.  Also  the  results  were 
included  obtained  by  the  project  authors  earlier  onboard  Prognoz,  Venera,  Vega  and 
Phobos  s/c. 

According  to  the  Project  2063  Work  Plan  the  work  has  been  performed  in  the  directions 
formulated  in  tasks  1-6. 

Task  1.  Selection  of  SEP  events  according  to  the  accepted  criteria  and  elaboration  of  the 
data  base  for  proposed  investigations. 

Task  2.  Detailed  examination  of  characteristic  decay  times,  spectra  and  anisotropy  of 
the  1  to  1 00  MeV  proton  intensities  in  SEP  events  observed  at  different  radial  distances 
over  azimuthal  range  available  for  observations. 

Task  3.  Looking  for  a  correlation  of  the  decay  time  of  particle  intensity  and  spectral 
index  with  the  speed  of  the  shock  associated  to  a  given  SEP-event. 
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Task  4.  Investigation  of  the  functional  form  (power  or  exponential  function)  of  particle 
intensity  during  the  event  decline. 

Task  5.  Studying  relationships  between  SEP  event  particle  fluences  and  both 
interplanetary  space  conditions  and  CME  parameters 

Task  6.  Analysis  of  variations  of  decay  times  over  solar  activity  cycle. 

All  principal  objectives  of  the  technical  schedule,  listed  in  the  Work  Plan,  were 
completed.  Sometimes  scope  of  activity  was  extended  according  to  the  logic  of 
investigations. 

According  to  the  Project  2063  Work  Plan  and  in  addition  to  it  the  work  has  been 
performed  at; 

•  selection  of  SEP  events  and  elaboration  of  the  initial  multispacecraft  data  base 
for  1977-1983  events; 

•  enlargement  of  the  data  base  due  to  IMP  8  1974-2001  data 

Statistical  analysis  of  the  next  features  of  decay  phase  in  SEP  events  has  been 
performed: 

•  Statistical  analysis  of  SEP  decays:  distributions  of  characteristic  decay  times, 
spectral  indices  and  total  durations  of  exponential  declines; 

•  Energy  spectra  of  SEPs; 

•  Periods  of  time  profiles  with  and  without  shocks; 

•  Functional  form  of  decays  of  SEP  events; 

•  Invariance  of  charged  particle  time  profiles  during  the  late  phase  of  SEP  events 
according  to  multi  spacecraft  observations. 

•  Correlation  between  decay  rate,  spectral  index  and  shock  speed 

•  Periods  of  quasistationary  conditions  in  interplanetary  space  according  to  the 
data  on  successive  SEP  events 

•  Influence  of  the  neutral  current  sheet  on  the  particle  fluxes  during  decay  phase 
of  SEP  events 

•  The  dependence  of  characteristic  decay  time  and  spectral  index  on  the 
flare  site 

•  Investigation  of  proton  fluences  and  adjacent  questions; 

•  Energy  dependence  of  characteristic  decay  time  in  exponential  decay  of  proton 
SEP  events 

•  Analisis  of  variations  of  decay  time  and  spectral  index  over  solar  activity  sycle 
The  results  of  investigatin  of  these  SEP  event  phenomena  are  described  below 
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Selection  of  SEP  events  and  creation  of  the  database 


One  of  the  main  stages  of  the  work  under  the  Project  was  creation  of  the  solar  energetic 
particle  (SEP)  events  for  most  prolonged  period  of  time.  It  was  possible  in  the  present 
investigation  on  the  basis  of  multi  spacecraft  measurements  to  embrace  a  period  from 
1974  to  2001,  that  is  more  than  two  solar  activity  cycles  (21-st  cycle  started  in  the 
middle  of  1976,  22-nd  -  in  the  middle  of  1986  and  23-th  -  in  the  middle  of  1996) 

For  systematic  analysis  of  the  decay  phase  of  SEP  events  at  first  selection  of  events 
according  to  the  shape  of  their  time  ptofiles  has  been  made:  only  the  events  with 
exponential  declines  were  took  into  consideration. 

Our  selection  criterion  for  SEP  events  has  only  a  limited  meaning.  It  is  clear  that  the 
intensity-time  profile  of  the  SEP  event  depends  not  only  on  the  number  of  accelerated 
particles,  but  is  also  influenced  by  four  factors.  These  are  (1)  the  angular  separation 
between  the  source  and  the  observers  magnetic  footpoint,  (2)  the  radial  distance,  (3)  the 
coronal  and  interplanetary  propagation  conditions  and  (4)  the  background.  Independent 
of  radial  distance  there  is  a  background  level  variation  of  a  factor  of  about  3  during  the 
solar  cycle,  beside  any  variation  of  background  due  to  a  preceding  event  (Kallenrode  et 
al.,  1992).  Using  the  adopted  criteria  the  database  of  gradual  and  impulsive  SEP-events 
related  to  flare/CME  associations  and  flares  only,  correspondingly,  was  done. 

Primary  data  base  consisted  of  122  gradual  and  100  impulsive  events  (see  Tables  lA  and 
2A  of  the  1-st  Year  Report).  In  the  Tables  there  are  the  lists  of  flares,  CMEs, 
interplanetary  shocks  and  spacecraft  data.  These  are  the  date  of  the  flare,  the  flare  onset 
and  location  in  Ha  and  the  Xt-burst  amlitude  from  Solar  Geophysical  Data.  CMEs  are 
presented  by  observed  locations  (position  angle  and  angular  width)  and  speeds.  These 
parameters  are  visible  projections  of  the  real  values  onto  the  sky  plane.  Under 
supposition  that  CME  originates  on  the  Sun  nearby  a  flare  site  (Reames  et  al.,  1996)  and 
moves  radially  from  the  Sun  we  reealculated  observed  speed  to  near  real  one  according 
to  procedure  discribed  in  (Daibog  et  al.,  1996).  This  value  is  also  shown  in  the  Table  lA. 
Shock  data  include  the  observation  points,  their  radial  distances  and  longitudes  relative 
to  the  flare  (azimuth)  and  average  (transit)  speeds. 

Unfortunately  observations  of  CMEs  began  only  from  March  1979  and  for  earlier  events 
we  could  use  only  data  on  interplanetary  shocks,  which  are  well  correlated  with  CMEs. 
The  spacecraft  data  in  the  database  include  locations  of  the  observation  points, 
characteristics  of  SW  and  IMF  at  them  and  the  coronal  footpoint  of  the  observer's 
magnetic  field  line. 

In  the  database  particle  data  are  the  values  of  fluxes  of  a  given  particle  specie  in  a  given 
energy  range  at  any  time  during  an  event  observation  at  a  given  s/c.  The  database 
contains  the  results  obtained  onboard  the  Helios  1  and  2,  IMPS,  ISEE3,  Prognoz, 
Venera,  SOHO  and  GRANAT. 

Patterns  of  time  profiles  of  protons  with  the  energies  >4,  >10,  >30  and  >60  MeV  are 
shown  in  Figs  la,b,c.  These  are  the  results  of  measurements  obtained  by  CPME 
instrument  onboard  IMP-8  satellite  in  1974-2001. 

The  Project  is  devoted  to  the  systematic  study  of  decay  phases  of  solar  energetic  particle 
(SEP)  events  in  the  inner  heliosphere  for  better  understanding  of  SEP  dynamics. 

For  systematic  study  of  decay  phases  of  solar  energetic  particle  (SEP)  events  all  flux 
rises  were  took  into  account,  independently  of  their  etymology:  whether  these  were 
flares,  corotating  or  ESP  events,  because  it  is  clear  that  in  space  weather  tasks, 
particularly  in  the  forecast  of  radiation  hazard  we  mu.st  not  limit  ourselves  to  «classical» 
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SEP  events,  associated  with  well  connected  flares.  That  is  why  the  primary  data  base 
was  enlarged  and  all  events  with  the  peak  flux  of  >4  MeV  protons,  exceeding  2-3  (cm' s 
sr)'  were  considered  independently  of  the  shape  of  time  profile,  decrease  rate  and 
anisotropy  of  particles  throughout  the  event. 

Only  those  exponential  parts  of  time  profiles  were  took  into  consideration  when  either 
exponential  dependence  continued  longer  than  12  hours,  or  intensity  decreased  more 
than  half  an  order  of  magnitude.  476  parts  of  all  time  profiles  with  exponential  declines 
were  selected. 


Fig.  la.  Pattern  of  proton  flux  time  profiles  for  energies  >4,  >10,  >30  and  >60  MeV 
beginning  from  Sept.  20,  1978.  Dashed  line  in  the  figure  top  shows 
selected  interval  of  the  exponential  decay  presented  in  the  Table  lA 
(Attachment). 
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Fig.  lb.  Pattern  of  proton  flux  time  profiles  for  energies  >4,  >10,  >30  and  >60  MeV  for 
sequence  of  events  beginning  from  Nov.3,  1997.  Dashed  lines  in  the  figure  top 
shows  selected  intervals  of  the  exponential  decay  presented  in  the  Table  1 A 
(Attachment). 
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Fig.  Ic.  Pattern  of  proton  flux  time  profiles  for  energies  >4,  >10,  >30  and  >60  MeV  for 
sequence  of  events  begirming  from  Apr.  9,  2001 .  Dashed  lines  in  the  figure  top 
shows  selected  intervals  of  the  exponential  decay  presented  in  the  Table  1 A 
(Attachment). 
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Statistical  investigation  of  the  decay  phase  of  SEP  events 

Over  the  whole  list  of  selected  events  the  values  of  characteristic  decay  times,  x,  of 
proton  Ep  >4  and  >10  MeV  fluxes  were  defined  at  the  exponential-shape  parts  of  time 
profiles.  These  values  are  presented  in  the  Table  1  of  Attachment.  This  table  includes  as 
well  the  values  of  peak  fluxes  of  protons  with  energies  >  10  MeV,  total  durations,  AT; 
spectral  exponent,  y,  for  the  beginning  of  exponential  decline  in  the  energy  range  4-60 
MeV  and  exponent  of  peak  flux  spectra  from  Catalogues  of  solar  proton  events  (Catalog 
of  Solar  Proton  Events,  ed.  by  Yu.I.  Logachev).On  the  basis  of  this  table  statistical 
investigation  of  the  values  of  characteristic  decay  times  of  all  exponential  parts  of 
declines  was  carried  out,  as  well  as  their  durations  and  spectral  indices.  Figs  2  a,b,c 
show  distributions  of  these  parameters. 

Accomplished  investigations  show  that  distributions  of  proton  flux  decay  rates  are  not 
strongly  dependent  on  their  energies,  at  least  in  the  energy  interval  4-10  MeV.  This  was 
analyzed  in  detail  when  we  considered  the  effect  of  temporal  and  spatial  invariance  in 
the  spectra  of  protons  and  electrons  in  the  major  SEP  events  (see  below).  By  virtue  of 
this  we  hope  in  the  future  to  consider  characteristics  of  the  decay  phase  only  for  one 
energy  and  chose  the  energy  4  MeV  as  the  most  representative  for  the  whole  set  of 
events.  However,  estimating  the  radiation  doze  through  the  whole  event,  it's  necessary  to 
take  into  account  different  energies  depending  on  the  thickness  of  protection  behind 
which  the  object  of  interest  is  placed. 
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Fig.2a.  Comparison  of  decay  time  distributions,  x,  of  >4  MeV  and  >10  MeV  protons. 
Both  are  normalized  to  100  events: 

1)  -  >4  MeV  protons,  2)  -  >10  MeV  protons. 

Fig. 2b  Distribution  of  proton  integral  spectrum  exponents,  y,  in  the  energy  range  4-60 
MeV  for  476  parts  of  exponential  declines  of  proton  fluxes. 

Fig.2c.  Distribution  of  SEP  events  on  the  duration  of  the  exponential  decay  of  the  >  4 
MeV  proton  flux 
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The  energy  spectrum  of  solar  cosmic  rays  is  the  most  important  characteristic  of  the 
acceleration  process  in  a  solar  flare  and  the  physical  conditions  in  the  acceleration 
region.  Various  acceleration  mechanisms  combined  with  different  conditions  in  the 
acceleration  region  give  different  energy  spectra  of  accelerated  particles.  In  particular, 
the  form  of  the  spectrum  and  the  total  amount  of  accelerated  particles  in  a  source  are 
closely  related  to  the  energetics  and  the  mechanism  of  the  flare  itself. 

It's  necessary  to  distinguish  between  current  instantaneous  spectra  during  an  event  and 
peak  flux  spectrum  that  to  some  extent  is  the  injection  spectrum.  We  examined  both  of 
them. 

The  maximal  values  of  proton  fluxes  of  different  energies  at  a  distance  of  1  AU  from  the 
Sun  were  used.  These  values  may  correspond  to  different  moments  of  time  for  particles 
of  different  energies.  If  the  emission  of  accelerated  protons  is  considered  to  be  pulse-like 
and  their  propagation  -  diffusive,  then  the  energy  spectrum  based  on  the  maximal 
observed  values  of  proton  fluxes  of  different  energies  will  be  the  spectrum  of  protons 
leaving  the  Sun  for  interplanetary  space,  i.e.,  be  emission  spectrum  of  solar  protons. 
Energetic  spectra  exponents  of  considering  events  are  limited  in  narrow  interval  with  the 
average  value  of  integral  spectra  exponent  1,7  ±0,8.  It  must  be  noted  that  energetic 
spectrum  of  protons  is  one  of  the  main  characterstics  defining  radiation  doze  behind 
different  protection  screens.  Most  dangerous  are  hard  spectra  with  small  values  of  y.  As 
all  events  that  have  occured  during  25  years  of  observation  (more  than  2  total  solar 
activity  cycles)  are  included  in  this  distribution,  one  can  definitely  say  that  no  serious 
deviations  of  energetic  spectra  from  their  average  values  can  happen.  In  the  present 
distribution  only  10%  of  spectra  have  values  of  y  less  than  1,  and  practically  all  of  them 
concern  the  low-power  events. 

Durations  of  the  exponential  parts  of  declines,  as  can  be  seen  from  the  Table  1,  are 
mostly  about  or  less  than  20  hrs.  Simultaneously  in  20%  of  events  durations  of  the 
exponential  parts  are  about  two  days  and  more.  As  a  rule  this  concerns  the  most 
powerful  events,  which  are  of  major  interest  both  from  the  point  of  view  of  particle 
propagation  in  interplanetary  medium  investigation  and  of  radiation  hazards  forecasting 
during  such  events. 

Periods  of  time  profiles  with  and  without  shocks 

One  of  the  Project  tasks  was  systematic  analysis  of  the  decay  phase  of  SEP  events 
associated  and  non-associated  with  the  shocks.  To  develope  understanding  of  processes 
of  propagation  and  losses  of  particles  and  a  role  of  interplanetary  shocks  in  forming 
intensity  time  profiles  during  the  decay  phase  at  different  dictances  from  the  Sun  and  at 
different  heliolongitudes  of  the  observer  multispasecraft  data  of  70-80-th  were  used,  as 
well  as  the  data  on  energetic  protons,  solar  wind  and  interplanetary  magnetic  field  from 
s/c  IMP-8. 

Interplanetary  shock  waves  and  CME-driven  shocks  strongly  influence  upon  flux  decay 
rates  and  on  SEP  spectra.  For  analysis  of  this  influence  distributions  of  t  and  y  values 
were  built  for  SEP  events  with  and  without  shocks.  These  distributions  are  presented  in 
the  Table  2  and  Figs.  3  and  4  (a,b,c). 

Figs.  3  and  4  (a,b,c)  and  the  table  permit  to  compare  distributions  of  the  values  of  t 
(Ep>4  MeV)  and  y  separately  for  time  intervals  accompanied  and  non-accompanied  by 
shocks.  One  can  see  that  shocks  steepen  declines  of  time  profiles  (x  decreases)  Average 
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values  of  t  during  shock  associated  and  non-associated  periods  are  10,9and  13.8  hrs. 
respectively. 

Distributions  of  y-s  are  more  compact  than  those  of  x,  there  is  no  considerable  deviations 
from  their  average  value  <y>  =  1,7.  Spectra  of  proton  fluxes  during  shock  associated 
periods  soften  (y  increases)  compared  with  non-shock  periods  (<y>  =  1,9  and  1,7. 
correspondingly).  Such  a  situation  show  that  acceleration  mechanisms  and  trapping  of 
particles  due  to  the  shock  is  more  effective  for  low  energy  particles. 


N 


Fig.  3.  Distribution  of  characteristic  decay  time,  x,  of  proton  >4  MeV  fluxes 

a) all  periods  with  exponential  declines 

b)  periods  without  shocks 

c)  periods  with  shocks 


Fig. 4.  Distribution  of  proton  integral  energy  spectrum  exponents,  y,  in  the  energy 
range  4-60  MeV.  The  value  of  exponent,  y,  relates  to  the  beginnings  of 
exponential  decline  periods 

a)  all  periods  with  exponential  declines 

b)  periods  without  shocks 

c)  periods  with  shocks 
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Special  analysis  of  total  durations  of  SEP  events  dependent  on  their  shock  association 
was  carried  out.  and  distributions  of  the  total  durations,  AT,  of  exponential  declines  with 
and  without  shocks  has  been  built.  Results  of  this  consideration  are  presented  in  Fig5. 
One  can  see  that  these  distributions  differ  insignificantly  for  AT>10  hrs,  while  for  AT 
<10  hrs  the  relative  number  of  shock  associated  events  is  30-  40  %  higher  than  without 
shocks.  This  can  be  due  to  the  fact  that  at  our  selection  criteria  short  duration  events 
often  represent  regions  of  shock-trapped  particles,  which  rapidly  pass  by  the  observer. 
Average  value  of  decays  without  shocks  AT  =  0,8  ±  02  days.  In  extreme  cases  duration 
was  longer  than  5-7  days.  These  events  have  a  large  value  of  characteristic  decay  time 
and  are  most  radiation  dangerous. 


Total  duration  of  decay,  hrs 


Fig.5.  Distributions  of  total  durations  of  decays 

-  accompanied  by  shock  (thin  line) 

-  non-accompanied  by  shock  (solid) 
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Functional  description  of  decays  of  SEP  events 


It  was  proposed  in  the  Project  to  investigate  the  functional  form  (power  or  exponential) 
of  the  SEP  event  declines.  This  is  of  a  great  importance  from  the  point  of  view  of 
adequacy  of  the  theoretical  models  considered. 

By  purely  diffusive  process  of  particle  propagation  flux  declines  must  be  power-law 
(J(t)=Af"  ,  where  n  =  3/2  for  classical  diffusion  in  the  three  dimentional  space),  by 
prevailing  of  process  of  convection  in  the  solar  wind  they  must  be  exponential-law 
(J=Joexp-t/T),  where  t  is  characteristic  decay  time  (Forman,  1970,  Jokipii,1972  ). 

As  was  shown  in  (  Daibog  et  al.,  2003a)  exponential-law  character  of  low  energy  (<  30 
MeV)  proton  flux  declines  are  observed  very  often  and  according  to  our  hypothesis 
could  be  a  criton  of  interplanetary  medium  state.  Really  a  rate  of  SEP  flux  decays  is 
defined  only  by  particle  spectrum  and  parameters  of  interplanetary  medium. 

As  was  told  above  decay  shape  is  of  exponential  form  if  convective  processes  prevail 
over  diffusion.  Then  i  could  be  described  as  x  =  3r/2Ur,  where  r  is  a  distance  of  the 
observer  from  the  Sun,  U  -  solar  wind  speed,  T  -particle  spectral  index  in  impulse 
representaion  J(p)  ~  p'^  (Forman,  1970;  Jokipii,  1972). 

In  the  energy  representation  J(E)~E‘'’',  that  gives  r=2y-l.  Forman  (1975)  and  Jokipii 
(1972)  have  obtained  for  F  the  next  expression:  r=2+a7,  where  a=(2MoC^+E)/(MoC^+E), 
E  -  particle  kinetic  energy,  MoC^  -  particle  rest  energy.  When  the  energy  value  is  small 
a~  2  and  F  =  2(l-i-y  ).  Somewhat  different  dependences  of  F  on  y,  obtained  in  these 
papers  do  not  influence  sufficiently  on  the  character  of  x  dependence  on  U  and  on  the 
spectral  index  y.  In  both  cases  x  is  reversely  proportional  both  to  U  and  y,  and  thus  for 
different  events  with  the  same  products  UF  the  values  of  x  will  be  the  same  as  well. 
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Fig.6.  Patterns  of  different  shapes  of  solar  proton  fluxes  registered  at  the  IMP  8 
s/c  after  their  peak  values: 

a)  the  event  with  the  exponential  decays  of  >10  and  >60  MeV  protons 
during  about  two  days 

b)  the  event  with  the  exponential  decay  of  low  energy  (>4  MeV)  protons 
and  a  power-law  for  the  high  energy  (>30  MeV)  protons. 

c)  the  event  with  a  long  (up  to  4  days)  exponential  decay  of  1  to  >10  MeV 
protons. 


In  fig.6  presented  are  the  patterns  of  events  with  different  functional  dependences  of 
proton  Huxes  on  time  during  their  declines.  Fig  6a  demonstrates  the  event  in  which 
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exponential  decay  of  protons  of  high  energy  (>10  and  >60  MeV)  was  observed  for 
nearly  two  days,  in  fig.6b  exponential-law  decline  is  valid  only  for  low  energy  protons 
(>  4  MeV)  while  for  high  energies  exponentiality  violates  and  for  >  30  MeV  protons  the 
functional  dependence  is  power-low.  Fig.6c  demonstrates  the  event  with  very  long 
exponential  decline  of  1  to  10  MeV  protons  for  the  event  of  June  7,  1979. 

Such  a  situation  (exponentiality  at  low  proton  energies  and  power-law  at  high  energies) 
is  in  agreement  with  theoretical  ideas:  exponential  decline  is  typical  under  predominance 
of  convection  over  diffusive  processes  Forman  (1975)  and  Jokipii  (1972),  which  is  more 
realistic  for  low  energy  particles,  though,  as  fig.6a  demonstrates,  sometimes 
exponential-law  decline  take  place  even  at  rather  high  energies  (60  MeV).  In  a  case  of 
event  presented  in  fig.6a  solar  wind  speed  was  600-700  km/s,  however  even  at  this  value 
of  solar  wind  speed  diffusive  processes  turned  to  be  supressed. 

It  was  shown  that  in  most  cases,  when  the  time  profile  following  the  maximum  is 
smooth  and  free  of  disturbances,  the  exponential  function  gives  the  best  approximation 
in  nearly  90%  of  the  cases  for  proton  energies  above  about  1  MeV.  Also  important  to 
note  is  that  this  exponential  decay  starts  very  soon  after  the  maximum,  not  only  in  the 
late  phase,  and  keeps  going  sometimes  for  one  to  two  weeks. 

The  comparison  of  experimental  values  of  characteristic  decay  times,  t,  with  those 
obtained  previously  in  different  theoretical  models  (Forman  1970,1971,  Owens,  1979) 
shows  that  theoretical  x  values  are  reasonably  close  to  the  fitted  slopes  in  nearly  50%  of 
all  cases  if  one  uses  the  average  solar  wind  speed  values, U,  measured  later  when  the 
corresponding  plasma  in  which  particles  were  convected  arrives  to  the  observer.  This 
might  be  a  surprisingly  high  number  bearing  in  mind  that  in'  rare  case  remains  U 
constant  for  a  sufficiently  long  time. 

For  mathematical  description  of  the  decay  phase  of  SEP  events  and  their  correlation 
with  IM  properties  it  was  supposed  that  the  declines  of  particle  fluxes  after  their  short- 
time  injection  are  defined  by  a  number  of  physical  processes,  the  main  of  which  are 
diffusion  (d)  and  convection  (c)  in  the  interplanetary  medium  and  adiabatic  deceleration 
(ad)  in  expanding  solar  wind.  Taking  into  account  these  three  factors  permits  to  define 
characteristic  decay  time  as 


1/X=  l/td+  1/Xc+  l/Xad 

Such  discrimination  of  processes  in  the  IM  permits  to  analyze  the  reasons  of  power  or 
exponential  declines  of  SEPs  after  their  maxima.  It  was  pointed  out  that  in  the  data  base 
created  for  this  Project  only  about  10  %  of  declines  are  described  by  the  power  law.  As  a 
rule  those  are  during  periods  of  weak  solar  wind  when  convection  may  be  neglected 
(1/xc  is  small  and  the  main  process  is  diffusion).  Such  a  picture  is  often  observed  for  high 
energy  particles  when  two  processes  -  convection  and  adiabatic  deceleration  may  be 
neglected,  especially  if  particle  energy  spectrum  is  hard. 

It  was  shown  that  the  power-law  decline  must  be  observed  also  in  a  case  of  adiabatic 
deceleration  of  particles  trapped  in  the  region  between  the  front  of  the  shock  going  away 
from  the  Sun  and  strong  magnetic  fields  near  the  Sun  (above  it  was  supposed  that 
particles  were  trapped  in  an  expanding  volume  element  in  the  solar  wind). 
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July  1985 


Fig.Vb.  A  pattern  of  the  power-law  flux  decline  in  July  9,1985  event  in  all  energy 
integral  channels  (>4,  >30,  >60  MeV). 

Upper  panel  shows  variations  of  solar  wind  speed. 


Sometimes  those  events  are  observed  in  which  low  energy  proton  (<  5-10  MeV) 
declines  may  be  described  by  the  exponential  law  (J  =  Jo  exp-t/r)  while  high  energy  (> 
10-30  MaB)  by  the  power  dependence  (J  =  Jo  t  The  patterns  of  all  three  kinds  of 
declines  for  high  and  low  energy  particle  are  presented  in  Figs. 7a  (exponential  declines 
for  all  energies),  7b  (power-law  declines  for  all  energies),  7c  (exponential  decline  for 
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low  energy,  power-law  decline  for  high  energy).  Upper  panels  represent  time  variations 
of  solar  wind  speed. 


November  1977 


Fig.7c.  A  pattern  of  the  22.1 1.77  event  with  the  exponential  decay  of  low  energy 

(>4  MeV)  protons  and  a  power-law  for  the  high  energy  (>30  MeV)  protons. 
Upper  panel  shows  variations  of  solar  wind  speed. 


It  must  be  noted  that  more  rapid  exponential  decay  is  more  favorable  from  the  point  of 
view  of  radiation  hazards  during  space  missions,  especially  for  technical  elements 
behind  thin  protection  screens:  solar  batteries,  optical  instruments,  detectors  of  scientific 
devices. 
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Invariance  of  charged  particle  time  profiles  during  the  late 
phase  of  SEP  events  according  to  multispacecraft  observations 

The  question  was  considered  about  existence  of  spectral  invariance  in  electron  events, 
similar  to  that  observed  previously  in  proton  SEP  events.  It  was  important  from  the  point 
of  view  of  searching  of  arguments  in  favor  of  electron  acceleration  by  CME  driven 
shocks. 

20  enhancements  of  4-13  and  13-27  MeV  protons  and  0.3-0.8  and  0.8-2.0  MeV 
electrons  were  selected  that  were  simultaneously  observed  at  Helios  1  and  Helios2  s/c  at 
the  distance  0.3- 1.0  AU  from  January  1976  to  March  1980  (detailed  description  of  the 
events  is  presented  in  (Daibog  et  ah,  2000)). 

The  analyses  showed  that  1 7  events  from  20  candidates  had  signs  of  invariant  spectra  of 
both  protons  and  electrons.  13  of  them  were  completely  invariant.  Some  characteristics 
of  these  events  are  placed  in  Table  1 . 

Table  1 


Parameters  of  SEP  events  with  invariant  spectra 


Date 

Ton 

X 

R,  AU 

AX 

Vsh,  km/c 

Heliosl 

Helios2 

Heliosl 

Helios2 

Heliosl 

Helios2 

19.09.77 

09.57 

W57 

0.62 

0.71 

E84 

E108 

760 

HA 

24.09.77 

05.50 

W120 

0.57 

0.64 

E25 

E48 

1160 

860 

22.11.77 

09.45 

W40 

0.66 

0.61 

W78 

W47 

470 

660 

01.01.78 

21.45  ■ 

E  06 

0.95 

0.93 

W34 

EOl 

1130 

950 

13.02.78 

01.16 

W20 

0.95 

0.95 

W65 

W69 

HA 

870 

08.04.78 

01.10 

W  11 

0.50 

0.51 

W49 

W20 

690 

720 

23.09.78 

09.44 

W50 

0.74 

0.72 

E71 

Elio 

800 

780 

11.12.78 

18.33 

E  14 

0.70 

0.75 

W41 

E05 

740 

16.02.79 

01.23 

E59 

0.98 

0.97 

E04 

E36 

HA 

950 

01.03.79 

09.55 

E58 

0.94 

0.93 

WIO 

E30 

980 

810 

09.03.79 

09.35 

E80 

0.91 

0.84 

Ell 

E52 

920 

590 

E  17 

0.46 

0,36 

imi 

E32 

880 

HA 

E36 

0.62 

0.75 

W34 

E18 

650 

900 

In  the  table:  the  date  of  the  flare  associated  with  the  SEP  event;  T  on  -  the  time  of  the 
flare  onset  in  Ha  (UT);  X.  -  is  flare  helioIongitude;R  -  is  the  observer's  heliocentric 
distance;  AX  -  is  the  observer's  longitude  relative  to  the  flare  (W  -  westward,  E  - 
eastward  from  the  observer's  point);  Vsh  -is  the  transit  shock  speed  at  the  observer's 
point  (DG  -  is  a  data  gap,  the  blank  -  is  absence  of  effect). 
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Two  examples  of  the  invariant  spectra  events  are  shown  in  Figs.  8a  and  8b. 


Fig  8a.  Time  profiles  of  proton  (a)  and  electron  (b)  intensities  in  the  event  of  Sept.  23 
(DOY  266),  1978.  J  -is  hourly  averaged  counting  rate. 

(a)  diamonds  -  4-13  MeV  protons,  crosses  -  13-27  MeV  protons  according  to 
Helios  1  data;  light  and  filled  squares  according  to  Helios  2  data, 
correspondingly. 

(b)  crosses  -  0.3-0. 8  MeV  electrons,  circles  -  0. 8-2.0  MeV  electrons  according 
to  Helios  1  data;  light  and  filled  triangles  according  to  Helios  2  data, 
correspondingly 

It  can  be  seen  that  late  after  the  time  of  the  intensity  maximum  and  the  passage  of  an 
interplanetary  shock  at  the  s/c  intensities  of  protons  and  electrons  with  the  same  energy 
are  cIo.se  in  magnitude  and  decrease  with  nearly  the  same  x  on  both  s/c.and  one  can  say 
about  "complete"  invariance  of  particle  spectra.  For  the  event  of  23.09.78  invariance  of 
proton  spectra  was  observed  over  a  longitude  interval  of  about  160°  (Reames  et  al., 
1996.  Reames  et  al.,1997).  Here  we  can  see  that  electron  intensity  profiles  have  the  same 
characteristics  as  the  proton  ones  and  this  continued  at  least  until  October  8, 
demonstrating  high  stability  of  conditions  in  a  very  large  spatial  region. 
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Fig.  8b.  The  same  as  fig.Sa  for  Feb.  16  (DOY  47),  1979  event 

The  event  of  16.02.79  was  accompanied  by  two  shocks.  The  first  one  is  associated  with 
16.02.79  flare  located  at  E59  at  1:44  UT,  the  second  -  with  the  flare  on  18.02.79  at  36 
UT  located  at  W13.  The  latter  was  impulsive  flare  and  caused  electron  enhancements  on 
both  s/c  .  The  decay  times  of  these  enhancements  are  considerably  less  than  late  in  the 
event  (5jl0  hrs  comparing  to  50  hrs  in  the  late  phase).  Thus  contribution  of  electrons 
from  this  additional  ejection  during  the  decline  of  the  whole  event  is  reduced  after 
midday  of  19.02.79. 

The  analyses  of  the  totality  of  events  shows  that  spectrum  invariance  late  after  maxima 
of  SEP  events  are  observed  for  both  protons  and  electrons.  Sometimes  x  for  electrons 
are  20-30%  larger  than  for  protons,  this  can  be  explained  by  better  storage  conditions  for 
electrons  compared  to  protons  in  the  region  between  the  shock  front  and  converging 
magnetic  field  near  the  Sun.  Invariant  spectra  seem  to  occur  in  all  gradual  SEP  events. 
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One  may  expect  that  the  more  to  the  west  from  the  line  Sun  center  -  flare  (i.  e.  from  the 
nose  of  the  shock)  the  s/c  is  located,  the  more  is  the  delay  of  invariance  onset  due  to 
waiting  of  the  arrival  of  field  lines  connected  to  the  eastern  flank  of  the  shock. 

In  a  case  of  events  coimected  to  the  eastern  flank  the  observer  just  from  the  onset  of 
event  locates  in  the  invariance  region  and  the  onset  of  this  regime  outstrips  the  shock 
arrival. 

It  was  found  in  Reames  et.al.(1996)  that  delay  (or  advance)  of  the  invariance  onset  times 
of  proton  intensities  increase  with  increasing  longitude  separation  of  spacecraft  from  the 
flare  location.  For  eastern  flares  this  is  explained  by  existing  of  time  intervals  between 
crossings  at  s/c  field  lines  connecting  to  the  quasi-parallel  region  on  the  eastern  flank  of 
the  shock  and  to  the  longitude  at  the  shock,  that  primary  arrival  at  the  observer,  while  in 
a  case  of  western  flares  s/c  is  connected  to  this  region  from  the  very  beginning. 


1  At .  ,ne 

i  ; 


Fig.  9.  Delay  (negative)  and  advance  (positive)  time  intervals,  At.  of  invariance  onset 
relative  to  the  shock  arrival  versus  corresponding  heliolongitude  distance, 

AX,  between  the  s/c  and  a  flare.  Straight  lines  join  the  points  corresponding  to 
simultaneous  measurement  at  Helios  1  and  Helios  2  s/c.  Figures  near  the 
points  -  transit  shock  speeds  at  the  observation  points 
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The  question  arises  whether  electrons  also  “feel”  this  difference  in  flare-spacecraft 
locations,  as  well  as  about  magnitudes  of  these  delays  and  advances.  We  analyzed  the 
proton  and  electron  profiles  and  were  convinced  that  in  13  of  17  events  with  invariance 
properties  it  was  possible  to  determine  the  time  of  onset  of  invariance  regime  both  for 
protons  and  electrons.  These  events  are  presented  in  the  Table  1 . 

In  Fig.9  plot  of  delay  and  advance  time  intervals  At  versus  corresponding  longitude 
distance,  ATi,  are  presented  for  all  events  considered  (Daibog  et.al.  2003b).  Here  left  side 
of  the  plot  concerns  the  events  from  the  eastern  relative  to  the  spacecraft  flares,  right  - 
from  the  western  ones. 

Positive  Ats  correspond  to  advances,  negative  -  to  delays.  In  a  case  of  protons  all  events 
associated  with  eastern  relative  to  s/c  flares  turned  to  be  in  the  quadrant  corresponding  to 
delays,  with  central  flares  -  become  invariant  together  with  the  arrival  of  the  shock,  most 
of  western  events  -  in  the  quadrant  with  advances. 

In  a  case  of  electrons  all  eastern  events,  but  one,  occupies  the  quadrant  with  delays,  all 
western  events  but  one  -  the  quadrant  with  advances.  One  can  see  the  tendency  of 
increasing  of  discussed  time  intervals  ,At,  with  increasing  AX. 

Additionally  for  illustration  of  this  conclusion  the  points  corresponding  to  observations 
of  the  same  event  at  different  points  are  joined  by  straight  lines.  The  values  of  shock 
transit  speeds  are  presented  near  each  point. 

Unfortunately  the  statistics  is  not  numerous,  measurements  were  made  at  different 
distances  from  the  Sun,  and  this  caused  large  scattering  of  the  points.  This  is  also  due  to 
considerable  variations  of  solar  wind  speed  and  possible  changing  of  angular  sizes  of 
shocks  at  different  distances  from  the  Sun. 

In  the  absence  of  diffusive  scattering  particle  intensities  in  the  region  of  invariant  spectra 
have  to  be  identical  all  the  way  from  the  Sun.  The  events  observed  by  Helios  1  and  2  at 
distances  0.3-0. 6  AU  and  possibly  near  the  Sun-Earth  line  where  considered.  It  was 
found  that,  in  spite  of  the  low  number  of  events  there  were  those  in  which  Helios 
intensities  became  nearly  equal  to  those  observed  at  IMP  8  after  the  shock  crossing  the 
respective  s/c.  (Daibog  et  al.  2001). Thus  the  effect  of  spatial  invariance  in  flare-shock 
associated  events  is  valid  not  only  in  wide  angular  regions  but  at  different  distances  from 
the  Sun  from  0.3  to  1  AU  as  well.  A  pattern  of  such  event  is  shown  in  Fig.  10.  That  was 
the  April  30,  1976  small  event  after  a  flare  at  S08W46,  20:47  UT:  relative  positions  of 
flare  and  the  s/c  were  similar  to  those  in  Sep  23,  1978  event.  Helios  1,  2  were  at  0.65 
AU  and  0.42  AU,  respectively,  and  at  longitudes  about  W160.  >4  MeV  proton 
intensities  at  all  three  s/c  become  very  close  at  the  time  after  shock  arrival  to  IMP-8. 
Marked  are  shocks  at  Helios  1  and  2  {IMP-8  shock  transit  speed  910  km/s).  The  angular 
distance  between  the  magnetic  foot  points  of  the  Helios  1  and  2  s/c  were  only  10  deg  but 
it  took  about  3  days  before  intensities  became  equal.  This  may  indicate  the  change  of 
particle  propagation  conditions  at  about  0.5  AU 
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Figure  10.  The  fluxes  of  proton  in  SEP  event  30  April  1976 


Correlation  between  decay  rate,  spectral  index  and  shock  speed 

Investigation  of  correlation  between  rates  of  decay  of  proton  >4  MeV  intensities,  which 
can  be  described  in  exponential  form  J  =  Jo  exp(-t/T),  spectral  index  in  integral 
representation  calculated  for  the  start  time  of  exponential  decline  and  associated  shock 
speeds  was  fulfilled  on  the  data  base  of  IMP-8  in  1974  -  2001 

A  plot  of  characteristic  decay  times,  t,  as  a  function  of  the  shock  speeds,  U,  has  been 
built  for  two  criterions  of  event  selection:  1  -  first  after  maximum  exponential  parts  of 
declines  (as  a  matter  of  fact,  this  includes  in  many  cases  shock  peaks);  2  -  parts  of 
declines  after  shock  peaks,  if  any.  The  first  criterion  concerns  the  time  just  after  shock 
crossing  the  observation  point,  the  second  one  could  be  attributed  to  particle  trapping, 
described  above.  The  table  of  events,  presented  in  the  Annual  report  was  used,  also  new 
table  2A  (see  attachment)  was  created,  containing  parameters  of  declines  for  second 
criterion  selection,  which  are  "mixtures"  of  the  full  table  data.  The  results  are  presented 
in  Figs  1  la  (including  shock  peaks)  and  1  Ic  (without  shock  peaks). 
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Fig.l  1.  Depence  of  characteristic  decay  time,  t,  and  spectral  index,  y  ,  on  the 
shock  speed,  U 

There  is  no  clear  functional  dependence  between  x  and  U.  (Fig.  1).  However,  in  spite  of 
very  large  scattering  of  points,  some  limiting  dependence  can  be  noted  which  means  that 
very  large  x  takes  place  only  for  small  values  of  U  and  that  high  shock  speed  is 
associated  with  small  x.  Figs  lib  (including  shock  peaks)  and  lid  (without  shock 
peaks)  show  correlation  between  spectral  index,  y,  and  a  shock  speed.  Again,  point 
scattering  is  very  large,  but  some  tendency  of  slight  increasing  y  with  the  shock  speed, 
U,  exists. 

The  task  was  to  try  to  look  for  a  correlation  of  the  decay  time  of  particle  intensity  with 
the  speed  of  the  shock  associated  to  a  given  SEP-event.  It's  essential  that  SEP-events 
were  observed  at  various  angular  distances  relative  to  the  shock  nose  and  that 
longitudinal  dependence  (as  well  as  radial)  of  shock  speed  must  be  taken  into 
consideration.  In  the  model  developed  by  Reames  et  al.  (1996)  particles  are  supposed  to 
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be  quasi-trapped  in  an  expanding  flux  tube  between  the  high  magnetic  field  near  the  Sun 
and  the  moving  high  magnetic  field  region  downstream  of  the  travelling  shock. 

Though  the  authors  suggest  declines  to  be  caused  only  by  the  expanding  volume  and 
adiabatic  deceleration,  as  a  matter  of  fact,  they  don't  take  into  account  reflection  from 
the  outwardly  propagating  shock  front.  If  this  would  be  considered,  we  may  hope  to 
obtain  a  correlation  between  shock  speed  and  decay  time  and  thus  to  confirm  the  above 
model. 

They  proceed  from  the  model  (Reames  et  al.  1 996)  that  describes  the  evolution  of  the 
particle  distribution  as  being  governed  by  transport  equation  in  mixed  coordinates 
(Ruffolo  1995)  and  write  down  a  full  Fokker-Plank  equation  taking  into  account 
streaming,  convection,  p-diffusion,  adiabatic  focusing  and  adiabatic  deceleraton. 

A  rigorous  derivation  of  the  decay  time  of  the  differential  intensity  would  require  a  full 
solution  of  transport  equation,  subject  to  appropriate  source  function,  initial  and 
boundary  conditions.  A  number  of  successive  simplifying  suggestions  (  ignoring  fresh 
particle  injection  from  the  receding  shock;  scatter-free  motion  in  the  whole  volume 
except  of  the  shell  mentioned  above;  negligible  cross  field  diffusion;  isotropic  particle 
distribution  with  negligible  gradient  along  the  flux  tube;  power  law  energy  spectrum 
with  constant  power  index  )  as  a  matter  of  fact  reduce  the  task  to  the  simple  isotropic 
diffusive  model  with  adiabatic  deceleration  in  the  expanding  solar  wind.  So  by  many 
suppositions,  simplifications  and  angular  averaging  there  was  obtained  the  solution  for  t- 
dependence  of  t,  that  is  intrinsic  for  any  diffusive  model.  So  in  the  framework  of  their 
suppositions  it  would  be  natural  to  expect  the  same  dependence  both  for  shock 
associated  and  non-associated  events,  because  adiabatic  deceleration  in  the  expanding 
solar  wind  takes  place  in  both  cases.  All  the  suppositions,  mentioned  above,  permit  to 
approximate  the  decay  of  the  differential  particle  intensity  by 

T  =  t/(2Y  +  2), 

where  t  is  the  time,  y  -  spectral  index. 

It  can  be  demonstrated  that  the  similar  result  may  be  obtained  just  formulating  the 
accepted  one  by  one  additional  suppositions,  simplifications  and  angular  averaging  as 
the  main  conditions  of  the  task  One  can  formulate  the  simplest  task  of  particle  motion  in 
an  expanding  one-dimensional  volume  and  assume  all  particles  moving  parallel  and 
reflecting  from  the  ends  going  away.  As  in  Reames  we  suppose  that  particles  at  any 
moment  are  uniformly  distributed  in  the  volume  and  their  momentum  spectrum  is 
expressible  in  the  form  p"^^.  Then  from  the  continuity  equation 

5F/at  +  a(F*ai/dt)/ai  +a(F*ap/at)/ap  =  o, 

taking  into  account  that  dp/p  =  -dl/l,  where  1  =  Iq  +Ut  is  the  volume  length,  U  is  the  shock 
speed,  we  obtain  for  x  =  -F^SF/dt,  at  t »  Iq/U 

X  =  t/(2y  -  2). 

As  this  consideration  is  absolutely  transparent  and  as  it  is  very  difficult  to  check  Reames 
et  al.  chain  of  simplifications,  one  can  suppose  as  a  right  this  solution.  As  a  matter  of 
fact,  both  of  them  very  approximately  concern  the  real  picture,  because  we  are  not  at 
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liberty  to  choose  the  radial  dependence  of  distribution  function,  which  has  to  come  from 
the  solution  of  a  complete  Fokker-Plank  equation.  But  the  essential  fact  is  that  in  both 
solutions  T  decreases  with  increasing  y.  Simultaneously,  as  mentioned  in  (Ruffolo,  1995) 
deceleration  alone  leads  to  an  intensity-decay  time  in  the  expanding  solar  wind 

T  =  C/(y-l), 

where  Xd  can  be  obtained  from 


l/xd  =  2Vsw/3r, 


Vsw  is  the  solar  wind  speed. 


One  can  see  that  y-dependence  of  this  formula  is  similar  to  the  above  ones,  but  x  does 
not  depend  on  t  and  inversely  proportional  to  Vsw 

Thus  a  question  of  correlation  between  x  and  y  is  closely  connected  with  the  task  of  x  - 
shock  speed  relation. 

Analyses  of  IMP  8  data  permitted  to  obtain  a  plot  of  decay  times  as  a  function  of 
spectral  indices  (Daibog  et  al.,  2003a).  The  parts  of  time  profiles  with  and  without 
shocks  were  considered  separately  and  the  results  are  presented  in  Figs  12a  (including 
shock  peaks),  12b  (without  shock  peaks)  and  12c  (without  shocks). 


Fig.  12.  Correlation  between  proton  >4  MeV  decay  time,  x,  and  spectral  index,  y 
a  -  for  first  criterion 
b  -  for  second  criterion 
c  -  without  shock 


The  tendency  of  decreasing  x  with  increasing  y  definitely  exists  in  shock  associated 
events  (especially  in  Fig.  12a).  But  some  limiting  dependence  for  the  parts  of  time 
profiles  without  shocks  must  be  noted  as  well.  Qualitatively  Figs  12b  and  12c  differ  for 
small  values  of  y  and  x.  However,  we  must  underline  that  this  situation  contradicts 
neither  to  pure  diffusive  model  nor  to  the  model  with  particle  trapping  between  the 
shock  front  and  the  Sun.  Most  probable  adequate  model  must  take  into  account 
continuous  particle  acceleration  by  shocks. 
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Periods  of  quasi-stationary  conditions  in  interplanetary  space 
according  to  sequences  of  SEP  events 


In  the  case  of  no  additional  particle  injection  during  the  decay  phase  of  event  decay  rate 
for  particles  of  definite  kind  and  energy  is  determined  by  the  state  of  the  interplanetary 
magnetic  field,  more  exactly,  by  its  inhomogeneity  spectrum  in  the  region  of  particle 
propagation.  Declines  of  low-energy  (1-10  MeV)  proton  fluxes  can  often  be  described 
by  exponential  law:  J  =  Jo  exp  (-  t/x),  where  t  is  the  characteristic  decay  time.  This  may 
be  explained  by  convective  effects  and  adiabatic  deceleration  of  particles  during  their 
propagation  in  interplanetary  medium. 

Statistical  investigation  of  proton  declines  indicated  (Daibog  et  al.,  2003c)  that  the 
characteristic  decay  times  in  SEP  events  change  in  a  wide  spread  from  1  to  100  hours 
with  the  average  value  of  x  «  14  hours.  Simultaneously  it  was  discovered  that  sometimes 
very  long  exponential  declines  are  observed  with  constant  x,  they  last  up  to  5  days  and 
longer.  Over  3  solar  cycles  (1971-2001)  about  one  hundred  declines  with  constant  x  for 
more  than  7  days  were  collected.  Even  longer  declines  are  observed  very  seldom, 
because  at  the  end  of  events  it  is  difficult  to  extract  solar  particles  from  the  background 
and  impossible  to  judge  about  the  state  of  the  interplanetary  medium  (IM)  in  the  absence 
of  an  indicator  of  this  state.  However,  sometimes  it  turns  out  that  after  some  SEP  event, 
the  second,  third  etc.  events  follow  one  by  one  and  thus  it  is  possible  to  check  the  state 
of  IM  through  a  longer  period  of  time. 

It  was  pointed  out  that  looking  at  successive  events  one  can  notice  the  existence  of 
rather  long  time  intervals  (two  weeks  to  one  month,  and  even  longer  in  a  few  cases) 
during  which  the  characteristic  decay  phases  of  events  in  these  periods  are  very  similar. 
Periods  with  short  decay  times  can  be  followed  by  long  ones.  Large  and  small  events  are 
observed  with  similar  declines  (Daibog  et  al.,  2003c). 

For  investigation  of  sequences  of  SEP  events  data  from  Prognoz  1,  2,  IMP  5,  8,  Helios 
1,  2,  VEGA,  GRAN  AT  and  SOHO  s/c  within  a  long  time  interval  from  1971  to  2001 
were  used.  Prognoz,  IMP,  and  GRANAT  are  Earth  orbiting  satellites  permanently 
located  at  1  AU,  SOHO  also  located  near  the  Earth  at  the  LI  libration  point  at  0,01  AU 
upstream  of  the  Earth,  while  Helios  and  VEGA  were  located  at  0.3-^1.0  AU  and  0. 7^1.0 
AU  correspondingly,  that  is  why  we  indicate  in  the  Table  2  the  radius  and  heliolongitude 
of  s/c  for  Helios  and  Vega  observations.  There  were  analyzed  exponential  decays  of 
proton  integral  >4  MeV  fluxes,  as  well  as  differential  1-5,  4,6-15,  1-20  MeV  fluxes  and 
selected  both  long  lasting  exponential  decays  and  sequences  of  events  with  similar  (with 
an  accuracy  20%)  values  of  x. 

Fig.  13  a,b  shows  two  patterns  of  such  a  decline.  Fig.  13a  shows  the  decay  of  fluxes  of 
proton  Ep  2-4,6  MeV  and  Ep  4,6-15  MeV  with  x  =  32±2  hours  through  two  weeks 
period.  Two  decays  of  >10  MeV  protons  in  November  1997  (Fig.  13b)  were 
characterised  by  x  =  15±1  hrs.  for  a  week.  One  of  the  outstanding  events  of  September 
23.  1978,  studied  earlier  by  Reames  et  al.  (1996,  1997)  and  Daibog  et  al.  (2000),  when, 
according  to  Helios  1  and  2  data,  an  exponential  decay  with  x  =  33  hrs  was  observed  for  about 
two  weeks. 
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according  to  Helios  1  and  2  data,  an  exponential  decay  with  x  =  33  hrs  was  observed  for  about 
two  weeks. 
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Fig.  13.  Two  patterns  of  double  SEP  events  (upper  panel:  the  events  of  August, 

15  and  25,  1980;  lower  panel:  the  events  of  November,  4  and  6,  1997)  with 
identical  decrease  rates  of  different  energy  proton  fluxes  as  observed 
at  IMP  8.  For  the  former  case  x  equals  32±2  hrs,  for  the  latter  -  x  =15-17  hrs. 


Figure  14  shows  a  pattern  of  the  sequence  of  3  events  of  the  total  duration  about  9  days. 
This  says  that  during  this  time  everywhere  in  the  inner  heliosphere,  at  least  nearby  the 
solar  equator  plane,  the  conditions  of  interplanetary  magnetic  Field  (IMF)  were  stable 
that  resulted  in  a  practically  constant  proton  decay  times  in  the  sequence.  It  turned  out 
that  such  sequences  were  observed  even  longer,  sometimes  up  to  one  month  and  more, 
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Near-Earth  Helfosphere  Data  (OMNI)  by  OMNIWeb 


Fig.  14.  Pattern  of  proton  flux  time  profiles  for  energies  >4  MeV  for  sequence  of  events 
beginning  from  Apr.  9,  2001. 
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Fig.  15.  Solar  proton  fluxes  registered  onboard  the  PROGNOZ  satellite  in  April-July 
1972.  Similar  declines  are  marked  by  dotted  line. 


According  to  Prognoz  1  and  2  data  from  April  to  August  1972,  that  is,  during  3.5 
months,  8  subsequent  events  were  observed  with  practically  the  same  declines  of  proton 
1-5  MeV  fluxes  (Fig.  15)  <t>  =  16±1  hours  (Vernov  et  al.,1972).  A  similar  situation  was 
registered  aboard  IMP  5  for  0.9-1. 5  MeV  protons  (Van  Hollebeke  et  al.,  1974).  This 
means  that  prior  to  two  extremely  powerful  flares  of  August  4  and  7,  1972,  the  Sun  did 
not  create  any  sufficient  disturbances  and  interplanetary  medium  for  more  than  three 
months  was  quasistationary  returning  to  this  state  after  all  flares  of  different  power  that 
happened  during  that  period. 


The  interpretation  of  this  fact  is  rather  obvious.  In  the  “old  diffusive  language”  we  might 
point  out  that  the  mean  free  path  was  constant  throughout  the  whole  period.  In  the  recent 
language  of  exclusive  shock  acceleration  one  would  formulate  this  as  a  similarity  of 
particle  holding  between  strong  magnetic  field  near  the  Sun  and  the  shock  front  and 
particle  acceleration  in  interplanetary  space. 

It  must  be  noted  that  due  to  solar  rotation  during  such  considerable  periods  of  time  a 
spacecraft  is  connected  by  magnetic  field  lines  to  different  regions  on  the  Sun:  both  with 
flare  sites  and  with  quiet  surface.  However,  the  interplanetary  medium  remains 
magnetically  homogeneous  in  a  wide  angular  range  that  manifests  itself  in  the  constancy 
of  X.  This  result  indicates  the  existence  of  long  intervals  comparable  with  the  period  of 
the  solar  rotation  of  quasi  constancy  of  conditions  that  determine  the  particle  lifetime  in 
the  inner  heliosphere.  If  long  duration  exponential  decay  is  observed,  it  means  that  the 
IM  state  is  stationary  in  a  wide  angular  interval  (one  day  observation  due  to  the  Sun 
rotation  equals  to  a  13.3°  longitude  interval).  One-week  duration  of  decline  with 
constant  x  corresponds  to  a  stability  of  conditions  in  about  a  quarter  of  inner  heliosphere. 

In  Table  2  exhibited  are  the  characteristics  of  reliably  identified  sequences  of  two  and 
more  events,  which  were  observed  as  a  rule  during  more  than  a  week  through  1971- 
2001,  that  is,  nearly  for  three  solar  activity  cycles.  One  can  see  from  the  Table  2  that  the 
values  of  x  range  from  5  to  50  hrs,  one  order  of  magnitude,  and  this  may  be  the  scale  of 
variation  of  IM  parameters  responsible  for  particle  propagation  (for  example  mean  free 
path  in  the  case  of  diffusive  propagation).  Average  value  of  x  along  all  sequences 
declines  equals  16  hours,  and  it  concerns  approximately  to  30%  of  the  whole)  time 
during  which  diagnostics  of  interplanetary  medium  with  the  help  of  sequences  of  events 
is  possible.  It's  necessary  to  note  that  this  value  of  x  coincides  with  <x>  =  14.7±3.0  hours 
obtained  in  analysis  the  whole  statistics  of  exponential  declines  during  1974-2001 
(Daibog  et  al.,  2003d). 

The  existence  of  stationary  conditions  in  IM  was  discussed  previously,  we  could 
mention  the  idea  of  “uniform  particle  reservoir”  by  Roelof  et  al.,  (1992)  and  McKibben 
et  al.,  (2001).  Also  it  was  discussed  in  connection  with  so  called  “superevents”,  which 
are  long  lasting  enhancements  of  the  interplanetary  particle  population  observed  initially 
in  the  inner  heliosphere  and  then  propagate  to  the  outer  heliosphere  (up  to  35  AU)  with 
speeds  of  about  800  km/s.  As  one  hypothesis  for  the  origin  and  storage  mechanism  for 
these  events  it  is  suggested  that  a  system  of  large  interplanetary  shock  waves  exists.  The 
collision  and  superposition  of  a  variety  of  flare  initiated  shocks  outside  of  Earth  orbit 
leads  to  the  efficient  particle  acceleration. 
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Table  2 


Onset  date  of 
series: 

year/mo/day 


1971.04.02 


1972.04.20 


1974.12.22 


1976.03.16 


1978.04.01 


1978.06.11 


1978.07.14 


1979.02.27 


1979.02.18 


1979.06.04 


1979.07.07 


1979.08.05 


1979.11.06 


1980.05.27 


1980.07.03 


1980.08.02 


1980.12.01 


1981.05.12 


1981.07.02 


1981.08.28 


1981.11.16 


1981.11.24 


1981.12.09 


1982.01.03 


1982.06.03 


1982.07.10 


1983.05.12 


1984.02.16 


1986.01.19 


1990.03.31 


1990.08.06 


1991.03.25 


1991.05.13 


1992.11.24 


1997.10.07 


1997.11.05 


1998.01.19 


1998.10.21 


1999.05.27 


2001.04.12 


S/c,  Ep  (MeV) 


IMP-8;  4-12,5 


IMP-8;  4-12,5 


H-l;4-  13 


IMP-8;  4-12,5 


IMP-8;  >4 


H-l,-2;4-  13 


IMP-8;  4-12,5 


H-l;4-  13 


IMP-8;  4-12,5 


lMP-8;  4-12,5 


IMP-8;  2-4,6 


H-l;4-  13 


lMP-8;  4,5-15,0 


IMP-8;  4,5-15,0 


IMP-8;  4,5-15,0 


H-l;4-13 


IMP-8;  4,5-15,0 


IMP-8;  4,5-15,0 


H-l;4-  13 


H-l;4-  13 


IMP-8,  >  4 


IMP-8,  >  4 


Vega-1;  4,5-13 


GRANAT;  1-20* 


GRANAT;  1-20* 


GRANAT;  1-20* 


GRANAT;  1-20* 


GRANAT;  1-20* 


SOHO;  0.7-6 


IMP-8;  4,6-15,0 


SOHO;  0.7-6 


SOHO;  0.7-6 


SOHO;  0.7-6 


IMP-8;  4,6-15,0 


hours 


11,2- 

-24,0 

14,6- 

-24,2 

11,0- 

-  14.0 

40,0- 

-42.0 

7.0- 

-  10,0 

22,0- 

-230 

23,0- 

-24,0 

13,7- 

-  14,5 

Heliocoordinates  of 
s/c  Helios  and 
Vega;  Distance 
from  Sun  (AU), 
heliolongitude 


Notes 


15,0-18,0 


17,5-22,0 


16,5  -  19,0 


R=0,9;  9=120“  HI 
R=0,9  9=170°  H2 


17,0- 

-  18,0 

11,5- 

- 13,0 

31,0- 

-37,5 

48,0- 

-53,0 

7,3- 

-  7,6 

13,5- 

-  16,0 

9,0- 

-  12,0 

18,0- 

-  19,3 

7,0- 

-  10,8  ± 

15,0- 

-  19,0 

9,0- 

- 13,0 

17,7- 

-  17,7 

13,2- 

-  15,6 

20,0- 

-24,5 

20,0- 

-23,0 

36,5- 

-38,0 

10,0- 

-  12,0 

11,8- 

-  12,5 

10.5- 

-  15 

15,0- 

-  18,0 

12.5- 

-  14.5 

11,7-14,5 


N  —  a  number  of  solar  events  in  the  sequence; 

AT  -  a  duration  of  the  sequence  from  the  start  of  the  first  to  the  end  of  the  last  day; 
mark  "S"  means  that  the  sequence  took  place  during  a  superevent 
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In  spite  of  conditional  conception  of  superevents  (see  for  example  (Muller-Mellin  et  al.. 
1986  and  Droge  et  al..  1992)),  the  most  important  moment  is  the  high  stability  of 
conditions  in  interplanetary  medium  in  the  inner  and  outer  heliosphere  for  long  periods 
of  time.  It  was  natural  to  suppose  that  individual  SEPs  happened  during  a  superevent 
have  some  common  featuteres  and  this  might  be  decay  time.  Droge  et  al.,  (1992) 
identified  16  superevents  between  1974  and  1986.  Those  cover  a  spectrum  of  sizes, 
ranging  from  large,  well-defined  events  (Muller-Mellin  et  al.,  1986)  to  less  obvious 
cases.  Really,  most  part  of  sequences  described  above  happened  during  periods  of  1 1 
superevents.  These  are  N  2,  4,  6,  7,  8,  9,  10,  11,  13,  15,  16  from  Droge  et  al,  (1992). 
They  are  marked  by  the  sign  "s"  in  Table  2.  The  September  23,  1978  event  with  two 
weeks  decay  also  took  place  during  a  superevent.  However,  other  sequences  in  1974- 
1986  do  not  belong  to  periods  of  superevents,  as  well  as  in  other  4  superevents  decay  times 
of  subsequent  SEPs  differed  from  event  to  event. 

Thus,  one  can  see  that  the  sequences  of  events  with  similar  t  occur  both  during  and 
between  superevent  periods.  This  means  that  situation  of  stable  propagation  or  trapping 
conditions  is  not  inherent  exclusively  to  the  periods  of  superevents.  However,  this 
investigation  is  important  for  space  weather  forecast.  If  a  SEP  event  has  occurred  during 
a  period  of  stable  high  particle  background  (an  order  of  magnitude  higher  than  quite 
value)  for  a  few  days,  then  for  the  following  particle  enhancement  one  may  expect  the 
decay  time  similar  to  that  in  a  previous  event  and  thus,  according  to  the  peak  intensity 
time,  the  total  event  duration.  This  also  permits  to  predict  the  total  radiation  dose 
throughout  the  event. 

Thus,  it  turns  out  periodically  that  either  the  interplanetary  medium  parameters 
determining  conditions  of  energetic  particle  propagation  (diffusion,  convection  and 
adiabatic  deceleration)  or  the  conditions  of  trapping  of  particles,  accelerated  in  flares  and 
at  coronal  and  interplanetary  shocks,  between  strong  magnetic  fields  near  the  Sun  and 
shock  front,  stay  quite  stable  for  long  periods  of  time.  The  investigation  of  such 
sequences  of  events  can  also  be  useful  for  forecasting  characteristics  of  following  SEP 
events. 


Influence  of  the  neutral  current  sheet  on  the  particle  fluxes 
during  decay  phase  of  SEP  events 


Crossing  of  the  neutral  current  sheet  during  observations  near  the  ecliptic  plane  is 
associated  with  a  recurrent  shock.  Speed  of  such  a  shock  is  about  400  km/s  and  its 
influence  on  the  structure  of  charged  particle  fluxes  may  be  significant. 

In  many  investigations  current  sheet  was  considered  as  an  obstacle  for  charged  particle 
propagation  as  this  is  a  region  where  co-rotating  standing  bow  shock  is  generated  which 
might  change  temporal  characteristics  of  particle  fluxes  at  both  sides  of  the  shock. 
Detailed  analysis  of  the  influence  of  current  sheet  on  the  features  of  SEPs  during  the 
rise  phase  was  performed  by  Kahler  et  al.(1996). 

Recent  observational  results  suggested  that  the  current  sheet  may  be  an  important  factor 
in  organizing  the  propagation  of  shocks  and  SEPs.  Kallenrode  (1993)  found  that  times  to 
SEP  onset  and  to  SEP  flux  maximum  were  more  delayed  for  E  ~  0.5  MeV  electrons  and 
E  ~  20  MeV  protons  when  the  observer  and  solar  source  region  were  on  opposite  sides 
of  the  current  sheet  than  when  they  were  on  the  same  side.  Shea  et  al.,  (1995)  reported 
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that  for  22  of  27  ground  level  events  (GLEs))  the  solar  source  regions  were  in  the  same 
magnetic  polarity  as  the  Earth.  This  was  taken  as  evidence  for  preferential  particle 
propagation  in  the  polarity  structure  of  the  solar  source  region.  Studies  of  interplanetar} 
disturbances  also  support  the  idea  of  an  important  role  for  current  sheets.  Henning  et  al. 
(1985)  found  solar  wind  speeds  to  be  slower  following  large  flares  in  the  same  source 
surface  polarity  as  the  Earth.  Since  their  speed  distributions  were  qualitatively  consistent 
with  the  results  of  Cane’s  (1988)  shock  studies,  their  work  also  suggests  that  shocks 
weaken  when  they  cross  a  neutral  sheet. 
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Days  of  September  1978 

Fig.  16a.  The  pattern  of  crossing  the  current  sheet  during  Sept.23,  1978  event. 
Time-dependences  of  proton  fluxes  (>  4  and  >10  MeV  -  two  lower 
panels),  as  well  as  solar  wind  speed  (upper  panel). 


Kahler  et  al.  (1996)  have  examined  134  large  SEP  events  to  look  for  the  effect  of  the 
heliospheric  current  sheet  on  the  propagation  of  SEPs  They  found  no  significant 
differences  in  the  onset  times,  risetimes,  or  peak  fluxes  between  SEP  events  with  same- 
polarity  source  regions  and  SEP  events  with  opposite-polarity  source  regions. 
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An  examination  of  the  rise  phases  of  several  individual  SEP  event  flux  profiles  als(' 
fails  to  show  the  modulation  expected  when  a  sector  boundary  is  crossed.  It  appears  that 
if  the  speed  or  strength  of  the  CME-driven  shock  depends  on  the  position  of  the  source 
surface  current  sheet,  that  modulation  does  not  significantly  affect  SEP  production,  at 
least  during  the  first  few  hours  of  shock  propagation  from  the  Sun. 


Days  of  August  1979 

Fig.  16b.  The  same  as  Fig.  12  for  Aug.  19,  1979  event  under  multiple  change 
of  magnetic  field  polarity. 

Influence  of  the  current  sheet  on  the  decay  rate  of  particle  fluxes  may  be  examined 
during  both  isolated  events  and  sequences  of  events  with  the  similar  declines. 

The  existence  of  sequences  with  coincident  decay  time  suggest  that  during  the  time 
interval  of  the  sequence  IM  was  in  a  steady  state  near  the  point  of  obervation.  However, 
as  this  point  (the  location  of  the  s/c)  travels  13.3°  a  day  relative  to  the  IMF  due  to  solar 
rotation,  the  observed  sequence  extends  to  a  considerable  angle:  about  90°  for  a  week- 
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long  series  and  1 80°  for  2-weeks.  Then,  within  this  region  the  propagating  conditions  for 
energetic  particles  will  be  similar  (for  exponential  decays,  the  values  of  x  vary  very 
slightly),  and  consequently,  the  IMF  assumes  invariant  characteristics,  that  is,  uniform 
over  heliolongitude  and  constant  with  time. 


14  19  24  29  3 


Days  of  Febrary  1984 

Fig.  1 6c.  The  pattern  of  crossing  the  sector  boundaries  during  two  subsequent 
events  in  Feb.  1984,  which  shows  no  considerable  influence  on  the 
fluxes  of  both  low  (>4  MeV)  and  high  (>30  MeV)  energy  protons 
(middle  and  lower  panels).  Upper  panel  represents  the  direction  of 
the  radial  component  of  the  interplanetary  magnetic  field. 

Does  this  stability  propagate  in  heliographic  latitude,  to  the  different  sides  of  the  neutral 
current  sheet?  To  answer  this  question  the  position  of  the  spaeecraft  (i.e  the  Earth)  with 
respect  to  the  neutral  sheet  was  studied  during  the  sequence  observed.  It  turned  out  that 
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many,  especially  long  sequences  begin  at  the  same  sign  of  the  IMF  (with  Earth  located 
above  -  or  below  -  the  neutral  sheet),  and  end  with  an  opposite  magnetic  field  polarity 
(the  Earth  crossed  the  neutral  sheet  in  between).  During  this  time  x  did  not  change.  An 
example  of  such  event,  where  the  crossing  of  the  neutral  sheet  takes  place  in  the  middle 
of  the  decay  phase,  is  the  event  of  1978.09.23,  and  the  value  of  x  did  not  react  to  the 
sheet  crossing,  only  the  particle  fluxes  exhibited  some  fluctuations  (Fig.  16a).  This 
observation  suggests  that  both  hemispheres  (above  and  below  the  current  sheet)  were 
identical  during  the  event  series  from  the  viewpoint  of  particle  propagation.  In  Fig.  16b 
another  pattern  of  event,  Aug.20,  1979  where  the  crossing  of  the  neutral  sheet  takes 
place  in  the  middle  of  the  decay  phase.  In  Figs.  16a  and  16b  also  time  variations  of 
plasma  density  and  solar  wind  speed  is  presented  (upper  panel). 
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Days  of  April  200 1 

Fig.  16d.  The  same  as  Fig.  14  for  Ep>  2  MeV  and  Ep>  4  MeV  proton  fluxes  for 
the  sequence  of  events  during  April  2001 . 
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Figs.  ]6c  and  16d  present  the  patterns  of  crossing  the  current  sheet  during  the  sequence.^ 
of  events  in  Feb.  1984  and  April  2001.  Upper  panels  are  time  variations  of  radial 
component  of  interplanetary  magnetic  field.  One  can  see  that  sector  boundary  does  not 
considerably  influence  either  on  the  proton  flux  or  decay  rate  of  the  events  in  the 
sequence. 

It  mus  be  noted  once  more  that  earlier  papers  arrived  to  contradicting  conclusions.  The 
existence  of  the  neutral  sheet  (in  the  solar  corona)  between  the  flare  and  foot  point  of  the 
magnetic  line  connected  to  the  s/c  is  interpreted  as  an  obstacle  for  particle  propagation 
between  the  flare  site  and  the  s/c. 


Investigation  of  proton  fluences  and  adjacent  questions 

Particle  fluence  during  the  event,  equally  with  the  energy  spectrum,  is  the  main  radiation 
characteristics  of  the  event  because  for  taking  into  account  total  radiation  influence  on 
the  studied  object  it's  necessary  to  know  the  total  number  of  particles  with  definite 
energy  fallen  into  it. 

Particle  fluence  is  defined  as  the  total  number  of  particles  of  definite  energy  fallen  down 
1  cm^  through  the  event.  It's  very  important  to  predict  the  value  of  fluence  before  the  end 
of  the  event,  which  can  be  as  long  as  a  week  or  more  after  power  flares.  Such  a 
prediction  in  extraordinary  situations  gives  a  chance  to  undertake  preventive  actions 
against  radiative  defeat  of  living  organisms  and  elements  of  space  devices. 

Study  of  fluences  was  fulfilled  on  the  basis  of  126  major  events  during  a  period  1974  - 
2001  that  is  for  more  than  two  solar  activity  cycles.  This  includes  the  whole  21-st  and 
22-nd  cycles,  as  well  as  the  major  events  of  the  23-rd  cycle.  Only  those  events  were  took 
into  account  in  which  proton  >4  MeV  flux  exceeded  50  particles/cm^ssr  which 
corresponds  to  the  average  fluence  of  about  10  particles/cm  sr  and  is  equivalent  to  the 
absorbed  radiation  dose  of  about  10  rad.  Such  a  dose  is  not  sufficiently  dangerous  but 
for  most  power  fllares  fluence  can  mount  to  10'°  part/cm^  that  is  1000  times  higher  and 
this  threaten  with  serious  radiation  defeat,  especially  taking  into  account  that  all  the 
energy  releases  in  a  thin  layer  equal  to  the  scatter  free  path  of  protons  with  definite 
energy  (  4  MeV  in  our  case),  which  corresponds  to  the  thickness  of  sensitive  layer  of 
solar  batteries.  Protons  of  higher  energy  penetrate  dipper  and  there  energy  releases  in 
larger  volume.  For  inner  parts  of  spacecraft  the  most  danger  are  protons  with  energies  > 
30-50  MeV. 

The  analyst  showed  that  approximate  prediction  of  proton  fluence  is  possible  according 
to  the  peak  flux  of  particles  that  can  be  described  by  the  simple  dependence 
lgF=algJmax.  This  is  well  seen  in  the  upper  panels  of  fig.  17a,b,c,d  where  fluence  is 
represented  as  a  function  of  proton  peak  flux,  Jmax,  for  the  set  of  energies:  4,  >10,  >30 
and  >60  MeV.  However  dependence  of  the  fluence  on  Jmax  somewhat  differ  from 
event  to  event  and  one  can  see  considerable  dispersion  of  points  in  the  figures.  This 
dispersion  is  explained  by  different  durations  of  events,  more  exactly  by  different  decay 
rates  of  particle  fluxes  after  maximums  of  events,  which  is  characterized  in  a  case  of 
exponential  declines  by  characteristic  time  x  defined  from  the  dependence  J(t)  =  Jmax 
exp(-t/T).  The  more  x,  the  slower  the  decay  and  the  more  the  fluence  value. 

For  taking  into  consideration  decay  rates  of  the  events  products  Jmax*x  were 
additionally  considered.  The  results  are  presented  in  fig.  17(light  circles)  and  one  can 
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see  that  scattering  of  points  is  considerably  reduced  comparing  with  the  filled  circles. 
Thus  this  parameter  permits  to  predict  the  value  of  fluence  over  the  whole  event  by  the 
first  1-2  days  from  the  beginning  of  declines.  This  is  most  important  for  large  gradual 
events  that  are  really  radiationally  dangerous. 


Fig.  17.  Dependence  of  proton  a)  for  Ep>4  MeV,  b)  Ep>10  MeV,  c)  Ep>30  MeV 

and  d)  >60  MeV  fluence  though  the  event  on  the  peak  intensity  Jmax(«)  and 
Jmax*i:  (o)  and  on  the  product  of  peak  intensity  and  characteristic  decay  time 

There  was  analysed  the  influence  of  some  parameters  of  interplanetary  medium  (average 
through  the  event  solar  wind  speed,  average  absolute  value  of  interplanetary  magnetic 
field  and  average  magnetic  field  decline  to  the  ecliptic  plane)  on  the  main  characteristics 
of  SEP  events;  fluence,  characteristic  decay  time  and  proton  peak  flux.  It  turned  that 
niegther  of  these  parameters  practically  influence  on  the  peak  flux  of  event,  while 
characteristic  decay  time,  x,  is  considerably  influenced  both  by  the  solar  wind  speed  and 
the  value  of  interplanetary  magnetic  field,  x  is  reversely-proportinal  to  both  of  these 
parameters  (decrease  rate  of  proton  intensity  increases  with  increasing  both  solar  wind 
speed  and  absolute  value  of  interplanetary  medium).  Patterns  of  these  dependenses  for 
protons  with  the  energy  >4  MeV  are  shown  in  Figs.  18  a,b.  Such  a  dependence  is 
somewhat  weaker  for  fluences,  however,  the  similar  conclusion  may  be  made  for  thier 
envelope. 

It  must  be  noted  absolute  absence  of  any  dependence  of  all  event  parameters  on  the 
magnetic  field  direction. 
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Fig.  18.  Dependence  of  the  characteristic  decay  time  on  average  through  the  event 
solar  wind  speed  (a),  average  absolute  value  of  magnetic  field,  I B I  (b). 


Energy  dependence  of  characteristic  decay  time  in  exponential 
decays  of  proton  SEP  events 

Decay  rates  of  particle  fluxes  in  different  events  are  practically  independent  from  the 
power  of  event.  Simultaneously  particles  of  different  energies  behave  differently  not 
only  during  the  flux  rise  until  its  maximum  but  also  during  decays  of  events.  As  a  rule  a 
value  of  X  decreases  with  increasing  of  energy,  though  a  small  number  of  events  show 
opposite  dependence. 

To  study  the  energy  dependence  from  a  large  number  of  exponential  declines  in  the  data 
base  (more  than  200)  the  events  with  exponential  decays  were  selected  for  which  x  could 
be  determined  within  an  accuracy  of  10%  for  high-energy  differential  channels  (4.6-15, 
15-25  and  25-48  MeV),  i.e.  sufficiently  powerful  events,  because  in  weak  events  protons 
with  the  energies  of  about  25  MeV  either  are  absent  or  their  fluxes  are  so  small  that  it's 
impossible  to  determine  characteristic  decay  time,  x,  due  to  comparable  event  and 
background  intensities.  Also  selection  criteria  in  this  investigation  was  that  the  time 
profiles  following  the  maxima  were  smooth  and  free  of  considerable  disturbances.  Here 
the  fluxes  of  4.6-15  MeV  proton  usually  reach  Ipart  /(cm^  s  sr  MeV)  or  slightly  less. 

109  such  events  were  found  during  the  period  1973-2001  (Daibog  et.  al.,  2003e),  and  the 
exponent  assuming  power-law  energy  dependence  x=CE““was  determined  (E  is  the 
proton  kinetic  energy).  Due  to  the  wide  energy  intervals,  the  errors  in  determining  x  was 
about  10%  due  to  not  strict  exponentiality  and  small  flux  variations.  The  accuracy  of 
determination  of  a  is  about  30%.  Presented  in  Fig.  19,  the  distribution  of  a  values  can  be 
considered  to  consist  of  three  separate  groups:  (1)  x  independent  of  proton  energy 
(-0.1<  a  <0.1,  40  events);  (2)  x  decreasing  with  energy  (a>  0.1,  57  events);  and  (3)  a 
small  group  with  x  increasing  with  energy  (a  <  -0.1,  12  events).  Thus  predominantly,  x 
is  either  independent  of  or  decreases  with  proton  energy.  A  constant  x  results  in  spatial 
and  temporal  invariance  in  the  spectra  of  energetic  particles  in  gradual  solar  events 
discussed  in  (Reames  et  al.,  1996;  Daibog  et.  al.,  2000))  when  nearly  identical  spectra 
are  seen  over  large  heliolongitude  region  some  time  before  or  after  shock  arrival 
depending  on  the  relative  positions  of  the  flare  and  observation  point.  We  assume  that  it 
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is  possible  to  interpret  the  second  group  (decreasing  t)  as  having  a  “memory”  of 
diffusion  propagation  at  the  early  stage  near  the  Sun.  Indeed,  diffusional  decay  time 
profile  (power-law  shaped)  is  proportional  to  where  (A,  is  scattering  mean 

free  path,  u  is  the  particle  speed)  is  diffusion  eoefficient  increasing  with  the  particle 
energy  under  reasonable  assumptions  about  A,  dependence  on  particle  energy.  Thus  even 
exponential  decays,  which  are  most  probably  related  to  convection  and  adiabatic 
deceleration  rather  than  to  diffusion, 


Fig.  19  The  distribution  of  the  values  of  a  for  109  SEP  events 

are  in  some  cases  evidently  are  influenced  by  diffusion  (Forman,  (1970);  Owens,  (1979). 
The  most  surprising  result  is  the  presence  of  the  group  with  negative  a.  This  group 
includes  events  with  practically  all  possible  values  of  x  from  5  to  30  hours.  Since  the 
background  was  subtracted,  this  means  that  a  negative  a  is  not  simply  an  instrumental 
effect  and  might  indicate  that  some  additional  exotic  effects  of  propagation  and/or 
acceleration  can  play  a  role. 


The  dependence  of  characteristic  decay  time  and  spectral  index 
on  the  flare  site 

Under  observation  of  solar  cosmic  ray  fluxes  near  the  Earth  the  observer  is  connected  to 
the  Sun  at  heliolongitude  about  W60  (it  can  vary  according  to  the  value  of  solar  wind 
speed).  If  the  flare  occurred  at  this  heliolongitude,  then  particles  generated  during  the 
flare  would  come  to  the  observer  by  the  shortest  way.  If  the  flare  took  place  at  other 
solar  meridian  it  would  be  possible  to  expect  another  rise  time  and  slower  decay  after 
peak  intensity. 

McCracken  et  al.  (1971)  studied  properties  of  decay  phase  of  few  solar  flare  events 
according  to  Pioneer  6-9  observations.  Following  them  we  write  the  cosmic  ray  density 
as  U(r,v|/,  T,  t).  In  this  function,  \\i  is  the  heliolongitude  of  the  intersection  of  the  nominal 
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Archimedes  spiral  through  the  point  of  observation  with  the  solar  surface,  y  is  reckoned 
westward  from  a  fixed  point  on  the  rotating  Sun,  where  r,  T,  and  t  are  the  distance  from 
the  Sun,  the  particle  kinetic  energy,  and  time  respectively.  Since  the  cosmic  radiation 
remains  associated  with  the  same  magnetic  tube  of  force  as  it  is  convected  out  of  the 
solar  system,  a  cosmic  ray  population  "co-rotates"  with  the  Sun.  That  is  the  population 
remains  associated  with  the  same  values  of  v|;  throughout  its  lifetime  in  the  solar  system. 
The  rate  of  change  of  U  with  time,  as  observed  by  a  spacecraft  is  given  by: 

dU/dt  =  5U/5v|/  di(//dt  +  5U/5t 

d\}ld\.  includes  the  effect  of  adiabatic  deceleration.  Clearly,  the  two  terms  in  this  formula 
will  add,  or  partially  cancel  depending  on  whether  the  observer  is  on  the  Eastern,  or 
Western  side  of  a  cosmic  ray  population.  Any  discussion  of  the  decay  phase  of  observed 
SEP  event,  therefore  requires  knowledge  of  the  dependence  of  U  upon  heliolongitude. 

As  was  discussed  above  exponential  declines  testifies  that  convective  processes 
dominate  over  diffusive  and  Forman  (1970)  has  shown  that  the  particle  density  will  vary 
in  time  as: 


U  =  Uoexp(-  2V(2  +  aY)t/3r) 
where  the  dependence  of  U  on  v|;  has  been  ignored. 

Then  taking  into  account  this  dependence  we  can  write  the  observed  variation  as 

dU/dt  =  aU/SiK  d\i//dt  -U/x,  x  =  3r/2V(2  +  ay) 

Approximating  the  total  change  in  U  to  an  exponential  in  which  dU/dt  =  -U/T,  where  T 
is  the  observed  time  constant,  then 

-1/T  =  U''au/av|;  dydt -1/x 

and  if  over  a  limited  range  of  \\t,  Lf'aU/ay  =  l/yo ,  then 
1/T=  l/\|;o  di(;/dt  +  1/x 

For  the  events  considered  by  McCracken  et.al.,  (1971)  dv(;/dt  =  0.54°hr  v|/o  is  about 
±30°  . 

Indirect  indication  for  such  a  dependence  was  also  obtained  in  (Dalla,  2003)  where  AT 
was  considered  as  a  function  of  i|;.  The  Helios  1,2  data  show  a  possible  trend  for  the 
total  duration  of  a  particle  event,  AT,  to  decrease  as  the  location  of  the  associated  flare 
changes  from  eastern  to  western  longitudes  with  respect  to  the  magnetic  footpoint  of  the 
detecting  spacecraft,  within  the  range  of  values  of  .A\[/  in  [-1 80,+90]. 

It  was  interesting  if  such  a  dependence  exists  for  the  whole  statistics  of  events 
considered  in  this  investigation. 

On  the  basis  of  more  than  200  events  the  dependence  of  the  characteristic  decay  time,  x, 
on  the  flare  heliolongitude  during  1974-1996  (two  complete  cycles  of  solar  activity)  was 
analyzed.  The  result  is  that  x  is  statistically  independent  from  the  flare  heliolongitude, 
while  the  particle  arrival  and  peak  time  in  some  cases  are  significantly  delayed  with  the 
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increase  of  the  distance  of  the  flare  site  from  the  optimal  heliolongitude  where  magnetic 
field  line  is  connected  to  the  observation  point.  The  dependence  of  t  on  the 
heliolongitude  is  presented  in  Fig.20  (upper  panel) 

McCracken  et  al.  (1971)  have  also  obtained  indications  for  the  events  coiisidered  that  at 
late  times  (>  2  days)  the  spectral  exponent,  y,  near  10  MeV  is  dependent  on  the 
heliolongitude  of  the  observer  relative  to  the  centroid  of  the  particle  population  injected 
by  the  flare.  This  effect  results  in  a  variation  in  spectral  exponent  over  the  range  2.0<  y  < 
4.5.  Again,  our  statistical  consideration  does  not  confirm  this  fact.  The  dependence  of  y 
on  heliolongitude  is  shown  in  the  Fig.20  ( lower  panel) 

This  says  that  after  SEP  event  maximum  particle  propagation  in  the  interplanetary 
medium  is  controlled  mainly  by  IM  conditions  in  the  neighborhood  of  observation  point, 
and  the  structure  of  coronal  magnetic  fields  does  not  influence  on  the  shape  of  time 
profile  during  this  phase. 


-90  -60  -30  0  30  60  90  120  150 
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Fig.20.  Scatter  plots  of  decay  times,  x  (upper)  and  spectral  index, 
y  (lower)  versus  flare  heliolongitude 

This  conclusion  is  also  important  from  the  prognostic  point  of  view  because  this 
facilitates  forecasting  of  the  total  radiation  dose. 
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Analysis  of  variations  of  decay  times  and  spectral  indices  over 
solar  activity  cycle 


Up  to  now  nobody  has  undertaken  systematic  consideration  of  such  the  question  and  it 
has  only  been  noted  as  a  tendency  of  more  rapid  intensity  decline  for  some  events  at  quiet 
Sun  periods. 


years 


years 


Fig.21a.  Temporal  variation  of  t  for  integral  channels  of  protons  (left  panel)  and 
differential  (right  panel) 


Figure  21  a  and  b  display  all  values  of  t  and  y  in  exponential  declines  according  to 
selection  criterion  all  along  the  time  from  1974  to  2001.  Here  i  is  characteristic  decay 


years  years 

Fig.21b.  Temporal  variation  of  y  for  integral(left  panel)  and  differential 
(circles;  1-4.6  MeV,  crosses:  4.6-25  MeV)  channels  (right  panel) 


time  of  proton  flux  at  energies  >4  MeV  (integral  channels)  and  0.5-48  MeV  (differential 
channels)  during  the  time  of  exponential  decline,  AT;  the  value  of  y  corresponds  to  the 
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initial  moment  of  the  exponential  decline  period.  The  variation  of  the  over  two  solar 
cycles  suggests  little  variation  with  energy,  although  in  individual  cases  they  can  differ 
considerably  which  means  that  y  can  change  during  the  time  AT. 

There  was  no  discrimination  according  to  the  character  of  event,  whether  it  was 
diffusive,  associated  with  shock  waves  or  other  phenomena.  In  Fig.  21a,b  all  possible 
kinds  of  particle  enhancements  are  represented  independently  from  their  nature,  and 
thus,  the  observed  spread  of  values  of  x  and  y  apparently  reflects  the  totality  of 
phenomena  in  the  interplanetary  medium  because  over  more  than  tree  solar  activity 
cycles,  most  probably,  everything  happened  that  could  happen  at  all.  Neither  x  nor  y 
does  experience  obvious  variation  with  the  level  of  solar  activity,  apart  from  the  minima, 
where  the  width  of  the  distributions  seems  to  drop,  possibly  reflecting  the  stationary 
character  of  the  interplanetary  medium,  to  which  it  returns  after  various  disturbances.  In 
solar  activity  minimum  there  are  less  disturbances  and  the  interplanetary  medium  can 
return  to  its  undisturbed  state.  At  maximum,  the  medium  is  disturbed  practically  all  the 
time.  The  number  of  declines  is  changing  from  cycle  to  cycle;  it  was  by  20%  higher 
during  the  21st  cycle  (1977-1987)  than  in  the  22"^*,  which  corresponds  to  the  relative 
power  of  these  cycles. 
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Conclusion 


Many  natural  features  of  behavior  of  solar  energetic  particle  fluxes  in  the  decay  phase  of 
SEP  events  were  ascertained  under  the  Project  investigation.  One  of  the  most  important 
practical  conclusions  is  that  it  is  possible  in  principle  to  forecast  the  total  number  of 
particles  through  the  event  (fluence)  on  the  base  of  the  beginning  of  exponential  flux 
decline.  This  may  be  useful  for  averting  breaking-down  elements  of  space  devices  and 
the  danger  of  radiation  defeat  of  space  crews  (forestalling  switching  off  the  equipment, 
delays  of  entrance  to  the  open  space,  change  for  radiation-safe  modules  of  spacecraft, 
etc.) 

Executed  investigations  are  perspective  from  the  point  of  view  of  their  further 
development.  A  number  of  revealed  dependences  can  not  be  explained  in  the  framework 
of  current  theoretical  models  and  need  an  adequate  theoretical  interpretation. 
Simultaneously  it  became  clear  during  the  work  under  the  Project  that  some  features 
must  be  considered  in  more  detail.  Thus  up  to  now  there  is  no  exhaustive  clearness  on 
the  complex  of  causes  influencing  on  the  characteristic  decay  time,  x,  which  changes 
from  event  to  event  by  one  order  of  magnitude  and  especially  interrelations  of  these 
causes.  It's  possible  that  available  experimental  data  are  not  enough  for  solution  and  new 
experiments  are  necessary,  which  would  be  founded  on  the  new  theoretical 
consideration  of  particle  acceleration  on  the  Sun  and  by  shock  waves,  as  well  as  their 
propagation  in  the  interplanetary  medium. 
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1.  Energetic  Electron  Spectra  in  Solar  Energetic  Particle  Events  Associated  and 
Non-Associated  to  Coronal  Mass  Ejections 

V.  Slolpovskii,  E.  Daibog,  G.  Erdos,  S.  Kahler ,  K.  Kecskemety  and  H.  Kunow 
Proceedings  of  27-th  International  Cosmic  Ray  Conference.  Hamburg.  Germany  2001 
V.  8,  p.  3454. 

We  considered  energetic  electron  spectra  in  solar  energetic  particle  (SEP)  events  after 
flare/coronal  mass  ejection  (CME)  associations  and  after  flares  only  using  measurements 
onboard  various  spacecraft  (s/c)  at  radial  distances  of  0.3-1  AU  during  the  21st  solar 
activity  cycle.  Statistics  included  about  a  hundred  events  of  both  types.  More  than  50 
events  of  each  type  were  observed  simultaneously  at  various  points  of  the  inner 
heliosphere.  Energy  spectra  in  the  range  0.3  to  3  MeV  were  generated  from  maximum 
flux  at  each  energy.  The  relationship  between  the  exponent  of  integral  electron 
spectrum,  y  and  CME  speed,  V  and  angular  width  was  considered  and  it  was  obtained 
that  electron  spectra  become  harder  with  increase  of  V.  The  best  correlation  between  y 
and  V  arrives  when  the  observer's  magnetic  footpoints  are  in  the  limits  of  the  CME 
angular  width.  In  this  case  the  best  fit  approximation  of  y(V)  looks  as  yocV'°  -‘'  with  y  and 
V  in  the  range  of  4.5-1 .5  and  500-2000  km/s,  correspondingly. 


2.  Decay  Phases  in  Gradual  and  Impusive  Solar  Energetic  Particle  Events 

E.  Daibog,  V.  Stolpovskii,  G.  Erdos,  S.  Kahler,  K. Kecskemety  and  H.  Kunow 
Proceedings  of  27-th  International  Cosmic  Ray  Conference,  Hamburg,  Germany  2001,  v. 
9,  p.  3631. 

We  have  studied  the  decay  phases  of  energetic  particle  (SEP)  events  associated  and  non- 
associated  to  coronal  mass  ejections  (CMEs)  and  interplanetary  (IP)  shocks  using 
multispacecraft  observations  in  the  end  of  1970s  and  1980s.  Here  several  examples  of 
successive  SEP  events  observed  simultaneously  at  widely  spaced  s/c  are  presented.  It  is 
pointed  out  that  often  during  long  time 

intervals  (up  to  one  month)  the  decays  of  >4  Me'V  proton  and  >0.3  MeV  electron 
intensities  are  very  similar  over  a  wide  range  of  angles  between  the  observer  and 
proposed  source  (up  to  60  o  and  more).  The  spatial  and  temporal  invariance  of  the 
particle  intensities  during  the  late  phase  of  large  proton  events  associated  to  the  fastest 
CMEs  and  CME-driven  shocks  is  valid  not  only  in  a  wide  angular 
regions  but  at  different  distances  from  the  Sun  from  0.3  to  1  AU. 


3 .  Invariance  of  charged  particles  profiles  during  the  late  phase  of  SEP  events 

according  to  multispacecraft  observations 

E.I.Daibog,  V.G.Stolpovskii,  S.W.Kahler 

“Kosmicheskie  issledovaniya”,  2003,  v.  41,  Jfo2,  p.  140-147 

Using  simultaneous  observations  from  Helios- 1  and  Helios-2  s/c  we  examined 
properties  of  spatial  and  temporal  invariance  of  E>4  MeV  proton  and  E>0.3  MeV 
electron  spectra  during  a  decay  phase  of  SEP  events  after  flare  -  CME  -  CME-driven 
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shock  associations.  Spectra  and  decay  times  of  eharged  particles  are  nearly  invariant 
over  large  angular  distance  from  the  flare  point.  According  to  the  west  or  east  position  of 
s/c  relative  to  the  shock  nose  we  observe  delay  or  advance  of  invariance  onset  times 
following  the  shock  arrival  times  increasing  with  angular  distance  growlh. 


4.  Statictical  characterictics  of  SEP  event  flux  declines  through  a  long  period  of 
time  (1974-2001) 

E.I.  Daibog,  S.  Kahler,  K.  Kecskemety,  Yu.I.  Logachev 
Izvestiya  RAN,  seriya  fizicheskaya,  v.  67,  Ns  4,  p.  482,  2003. 

On  the  basis  of  statistically  homogeneous  data  from  IMP  8  statistical  analysis  of 
characteristic  decay  times,  x,  of  proton  >  4  MeV  fluxes  was  performed  during  time 
intervals  with  the  exponential  flux  declines.  For  these  intervals  the  exponents  of  integral 
energy  spectra,  y,  in  the  energy  range  4-60  MeV  were  determined  as  well.  Correlation  of 
T  and  Y  with  interplanetary  space  conditions  under  presence  or  absence  of  shock  w'aves 
were  analyzed. 

5.  Periods  of  Quasi-Stationary  Conditions  in  Interplanetary  Space  According  to 
Sequences  of  SEP  Events 

E.I.  Daibog,  Yu.I.  Logachev,  S.  Kahler,  K.  Kecskemety 
COSPAR02-A-00646  abstract ,  Houston,  USA,  2002 

Characteristic  decay  times  of  particle  fluxes  in  SEP  events  as  a  rule  vary  considerably 
from  event  to  event.  We  discovered,  however,  that  sometimes  the  successions  of  events 
with  very  close  decay  times  were  observed  during  long  time  intervals  (up  to  month,  and 
longer  in  few  cases).  The  values  of  characteristic  decay  times  differed  in  different 
successions. 

Such  successions  of  events  took  place  during  practically  all  superevents  in  the  21-st 
cycle  of  solar  activity.  However,  we  can  see  few  time  intervals  characterized  by  similar 
flux  declines  between  periods  of  superevents  as  well.  Thus,  we  can  state  that  it  turns  out 
periodically  that  either  the  interplanetary  medium  parameters  determining  conditions  of 
energetic  particle  propagation  (diffusion,  convection  and  adiabatic  deceleration)  or  the 
conditions  of  trapping  of  particles,  accelerated  in  flares  and  at  interplanetary  shocks, 
between  strong  magnetic  field  near  the  Sun  and  shock  front,  are  rather  stable  during  long 
periods  of  time. 

Investigation  of  such  successions  of  events  may  be  useful  from  prognostic  point  of  view. 


6.  Statistical  Study  of  Spectral  Characteristics  of  SEP  Event  Flux  Declines  during 
1974-2001. 

E.I.  Daibog,  Yu.I.  Logachev,  S.  Kahler,  K.  Kecskemety,  S.  McKenna-Lawlor 
COSPAR02-A-00546  abstract ,  Houston,  USA,  2002 

Interplanetary  medium  is  not  an  empty  space  for  propagating  of  previously  accelerated 
particles,  but  it  itself  influence  on  particles  changing  their  energy  (accelerating,  as  well 
as  decelerating  them)  due  to  interaction  with  magnetic  field  inhomogeneouties.  These 
processes  manifestate  in  energy  spectra  of  SEP  events.  Here  integral  energy  spectra  of 
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protons  in  the  energy  range  4-30  MeV  are  investigated  on  the  basis  of  a  homogeneous 
data  set  obtained  by  the  CPME  instrument  aboard  the  IMP  8  satellite  between  1974  and 
2001  for  all  SEP  events  with  >4  MeV  proton  flux,  exceeding  2-3  part/cm  s  sr. 

Another  important  characteristic  of  SEP  events  is  the  rate  of  decrease  of  particle  flux, 
which,  as  well  as  peak  flux  time,  is  an  integral  character  of  the  interplanetary  medium 
within  a  considerable  volume,  surrounding  an  observation  point. 

It  is  shown  that  the  average  values  of  spectral  index,  y,  and  characteristic  decay  time,  x, 
are  changing  along  with  the  solar  activity  phase. 

The  distributions  of  y  and  x  are  obtained  in  SEPs  with  and  without  shocks  and  during 
different  phases  of  events:  just  after  peak  flux  and  late  after  maximum. 


7.  Sequences  of  SEP  Events  as  an  Instrument  of  Indicating  Quasi-stationary  States 
of  Interplanetary  Medium 

Daibog  E.I.,  Logachev  Yu.I.,  Kahler  S.W.,  Kecskemety  K. 

“Kosmicheskie  issledovaniya”,  2003,  in  press 

Time  profiles  of  energetic  solar  particles  observed  in  the  interplanetary  medium  (IM) 
after  solar  flares  is  formed  to  a  great  extent  by  the  structure  of  interplanetary  magnetic 
field  (IMF)  and  its  inhomogeneouties  moving  from  the  Sun  with  solar  wind  speed. 
Decay  of  particle  flux  after  its  maximum  as  a  rule  is  independent  from  the  flare  power 
and  formed  mainly  by  the  propagation  conditions  in  the  IM.  Under  purely  diffusive 
propagation  of  solar  particles  flux  decline  is  described  by  the  power  law.  However  very 
often  it  turns  out  exponential-law  that  indicated  considerable  influence  of  convection 
and  adiabatic  deceleration  during  particle  propagation  in  the  expanding  solar  wind. 

In  the  present  paper  SEP  events  with  long  duration  exponential  declines  are  considered 
as  well  as  sequences  of  events  with  similar  exponential  declines  lasting  one-two  weeks 
and  longer.  The  existence  of  sequences  with  coincident  decay  time  suggests  that  during 
the  time  interval  of  the  sequence  the  IM  was  in  a  steady  state  near  the  point  of 
observation  and  due  to  solar  rotation  permits  to  judge  about  the  heliolongitudinal 
extension  of  the  homogeneous  region  in  the  IM.  As  a  rule  this  extension  over 
heliolongitude  is  from  a  quarter  to  a  half  of  the  near  ecliptic  space  and  even  occupies  the 
whole  near  the  Sun  space  in  few  cases  of  sequences  of  similar  proton  declines  lasted  for 
about  a  month  and  longer. 


8.  Some  Statistical  Properties  of  the  Decay  Phase  of  SEP-events 

K.  Kecskemety,  E.I.  Daibog,  S.  Kahler,  Yu.I.  Logachev 

Proceedings  of  28-th  International  Cosmic  Ray  Conference,  Tsukuba,  Japan,  Jul-Aug, 
2003,  in  press 

Generalized  parameters  characterizing  the  state  of  the  interplanetary  medium  (IM) 
include  the  functional  form  and  the  rate  of  decline  of  charged  particle  flux  in  solar 
energetic  particle  (SEP)  events.  The  shape  of  the  particle  flux  decline  is  of  particular 
importance:  power-law  time  dependence  indicates  the  dominance  ot  diffusive 
propagation,  whereas  exponential-law  decline  emphasizes  convection  transport. 
Depending  on  the  solar  wind  speed  both  exponential  (at  <10  MeV)  and  power-law  (>30 
MeV)  declines  can  be  present  in  the  same  event.  A  statistical  investigation  of  SEP 
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events  extended  for  a  long  period  suggests  that  about  90  %  of  SEP  decays  are 
characterised  by  exponential  declines.  Distributions  of  the  total  durations,  T,  of 
exponential  declines  with  and  without  shocks  differ  insignificantly  for  T>10  hrs.  Values 
of  the  decay  time  obtained  theoretically  are  reasonably  close  to  the  fitted  slopes  in  nearly 
half  of  all  cases  if  one  uses  the  average  V  solar  wind  speed  values  measured  when  the 
corresponding  plasma  in  which  particles  were  convected  arrives  to  the  observer. 
Dependences  as  the  variation  of  the  decay  time  with  energy  and  angular  distance  of  the 
observer  from  the  flare  heliolongitude  are  considered  as  well. 


9.  Statistical  characteristics  of  SEP  event  flux  declines  through  a  long  period  of 
time  (1974-2001) 

Daibog  E.I.,  Logachev  Yu.I.,  Kahler  S.W.,  Kecskemety  K. 

Abstracts  of  18-th  ECRS,  Moscow,  SH04P,  2002. 

The  rate  of  decrease  of  particle  flux,  as  well  as  peak  flux  time,  is  an  integral  character  of 
the  interplanetary  medium  within  a  considerable  volume,  surrounding  an  observation 
point.  We  performed  a  statistical  analysis  of  the  rates  of  decrease  in  all  SEP  events  with 
>4  MeV  proton  flux,  exceeding  2-3  part/cm^  s  sr,  on  the  basis  of  a  homogeneous  data  set 
obtained  by  the  CPME  instrument  aboard  the  IMP  8  satellite  between  1974  and  2001. 
We  consider  the  intervals  of  intensity  time  profiles,  during  which  declines  were 
characterized  by  exponential  law  typical  for  <50  MeV  protons  of  most  SEP  events.  The 
distribution  of  characteristic  decay  times,  x,  of  >4  MeV  and  >10  MeV  protons  for  355 
events  has  been  built. 

The  comparison  of  a  complete  set  of  events  from  21  and  22  solar  activity  cycles 
indicates  that  both  the  total  number  of  events  and  the  average  values  of  t  differed  in 
these  cycles.  Integral  energy  spectra  of  proton  fluxes  in  the  energy  range  4-30  MeV  were 
also  obtained  for  the  above  time  periods.  In  addition,  we  analyze  the  correlation  of  r,  the 
durations  of  exponential  declines  and  spectral  exponents  with  interplanetary  medium 
conditions,  in  particular  with  CMEs  and  interplanetary  shocks. 


10.  Invariance  of  proton  time  profiles  during  the  late  phase  of  SEP  events: 
multispasecraft  observations. 

Daibog  E.I.,  Stolpovskii  V.G.,  Kecskemety  K. 

Pisa  Solar  Wind  10  conference,  Pisa,  Italy,  2002,  abstract  SIV  43. 

Using  simultaneous  observations  from  Helios- 1,-2  and  IMP-8  s/c  we  examined  the 
properties  of  spatial  and  temporal  invariance  of  Ep>4  MeV  proton  and  Ee>0.3  MeV 
electron  spectra  during  the  decay  phase  of  SEP  events  after  flare-CME  -  CME-drived 
shock  associations.  The  spectra  and  decay  times  of  charged  particles  are  nearly  invariant 
over  large  angular  distances  from  the  flare  point.  According  to  the  west  or  east  position 
of  s/c  relative  to  shock  nose  we  observe  delay  or  advance  of  invariance  onset  times 
following  the  shock  arrival  times  increasing  with  angular  distance  growth. 
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11.  Diagnostic  of  the  state  of  the  interplanetar}'  medium  according  to  the  decay 
phase  of  SEP  events. 

Daibog  E.I.,  Kecskemety  K. ,  Kahler  S..  Logachev  Yu. I.  EGS-AGU-EUG  Joint 
Assembly  Nice,  France,  2003  Gephysical  Research  Abstracts,  Vol.  5,  02204,  2003 

Generalized  parameters  characterizing  the  interplanetary  medium  (IM)  state  are  the 
shape  and  the  rate  of  decline  of  charged  particle  flux  in  solar  energetic  particle  (SEP) 
events  (just  as  the  temperature  of  persons'  body  to  some  extent  is  a  generalized 
characteristics  of  his  state  of  health). 

The  shape  of  the  particle  flux  decline  is  of  particular  importance:  power-law  time 
dependence  says  about  prevalence  of  diffusive  propagation  of  particles,  exponential-law 
decline  about  convective  processes.  Depending  on  the  solar  wind  speed  both  kinds  of 
declines  may  be  realized  in  the  same  event:  exponential  for  low  energy  (<  10  MeV) 
particles  and  power-law  for  high  energy  (>30  MeV)  ones. 

Statistical  investigation  of  SEP  events  during  a  long  period  of  time  (from  1968  to  2001) 
shows  that  90  %  of  SEP  decays  are  characterised  by  exponential  declines. Distributions 
of  the  values  of  characteristic  decay  time,  x,  shows  that  most  part  of  the  time  during 
which  diagnostics  is  possible  the  interplanetary  medium  is  in  a  "basically  disturbed" 
state  characterised  by  t=16  ±  6  hrs.  This  permits  to  say  about  predominant  quasi  stable 
state  of  the  interplanetary  medium,  which  takes  about  one  half  of  the  total  time. 

As  SEP  events  usually  last  few  days,  the  rest  time  such  a  diagnostics  is  inaccessible. 
However,  often  SEPs  follow  one  by  one  and  all  events  in  the  sequences  have  similar 
decay  times,  x,  which  indicates  stability  of  conditions  in  the  IM  during  the  hole 
sequence.  This  permits  to  analize  the  IM  state  not  only  for  2-4  days  of  individual  event 
but  during  2-3  weeks  and  longer.  The  most  prolonged  sequence  and  IM  stable  state  was 
observed  in  1972  by  IMP-5  and  Prognoz-1.  Duration  of  this  period  took  above  2 
months,  i.e.  more  than  2  full  Solar  rotations. 


12.  Successions  of  Solar  Particle  Events  as  an  Instrument  for  Diagnostics  of 
Interplanetary  Medium  Conditions. 

E.I.Daibog,  Yu.I.Logachev 

Abstract  of  the  report  at  the  conference  of  CIS  and  Baltia  countries  "Actual  problems  of 
physics  of  solar  and  stellar  activity"  Nyzhnii  Novgorod,  June  2003,  p.  86 


The  values  of  characteristic  decay  times  of  particle  fluxes  in  SEP  events  as  a  rule  vary 
considerably  from  event  to  event.  We  discovered,  however,  that  sometimes  the 
successions  of  events  with  very  close  decay  times  were  observed  during  long  time 
intervals  (up  to  month,  and  longer  in  few  cases).  The  values  of  characteristic  decay 
times  differed  considerably  from  succession  to  succession. 

Such  successions  of  events  took  place  during  practically  all  superevents  in  the  21-st 
cycle  of  solar  activity.  However,  we  can  see  few  time  intervals  characterized  by  similar 
flux  declines  between  periods  of  superevents  as  well.  Thus,  we  can  state  that  it  turns  out 
periodically  that  either  the  interplanetary  medium  parameters  determining  conditions  of 
energetic  particle  propagation  (diffusion,  convection  and  adiabatic  deceleration)  or  the 
conditions  of  trapping  of  particles,  accelerated  in  flares  and  at  interplanetary  shocks, 
between  strong  magnetic  field  near  the  Sun  and  shock  front,  are  rather  stable  during 
long  periods  of  time. 
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Decay  Phase  of  Solar  Proton  Fluxes  on  the  Base  of  Statistical  Investigation  of 
About  300  SEP  Events 

E.I.Daibog,  Yu.I. Logachev 

Abstract  of  the  report  at  the  conference  of  CIS  and  Baltia  countries  "Actual  problems  of 
physics  of  solar  and  stellar  activity"  Nyzhnii  Novgorod,  June  2003,  p.  90 

The  interplanetary  space  is  not  a  passive  medium  for  propagating  previously  accelerated 
particles,  but  in  itself  influences  the  particles  by  changing  their  energies  due  to 
interactions  with  magnetic  field  inhomogeneities.  In  this  report  the  integral  energy 
fluxes  of  protons  in  the  energy  range  4-30  MeV  are  investigated  on  the  basis  of  two  data 
sets.  The  first  is  a  homogeneous  set  of  data  obtained  by  the  CPME  instrument  aboard  the 
IMP  8  satellite  between  1974  and  2001.  Selection  criterium  is  >10  MeV  proton  fluxes 
exceed  1 0  particle/cm^  s  sr.  The  second  consists  of  lower  energy  proton  fluxes  in  the 
range  of  0.75-6  MeV  measured  by  the  LION  instrument  aboard  SOHO  between  1995 
and  2001. 

Important  characteristics  of  SEP  events  include  the  rates  of  decrease  of  particle  flux, 
which,  as  well  as  peak  flux  time,  is  an  integral  characteristic  of  the  interplanetary 
medium  within  a  considerable  volume,  surrounding  an  observation  point.  It  is  shown 
that  the  average  values  of  energy  spectral  index,  y,  and  characteristic  decay  time,  x,  are 
all  changing  with  the  solar  activity  phase.  The  distributions  of  y  and  t  are  obtained  in 
SEPs  with  and  without  shocks. 


14.  Statistical  Properties  of  the  Decay  Phase  of  SEP-events 

Kecskemety  K.,  Daibog  E.L,  Kahler  S.W.,  Logachev  Yu.I. 

Abstracts  of  International  Solar  Cycle  Studies  Symposium  2003:  "Solar  variability  as  an 
input  to  the  Earth's  environment,  June23-28,  Tatranska  Lomnica,  Slovak  Republic, 
Abstract  VIII  7p,  p.l08 

Generalized  parameters  characterizing  the  state  of  the  interplanetary  medium  (IM) 
include  the  functional  form  and  the  rate  of  decline  of  charged  particle  flux  in  solar 
energetic  particle  (SEP)  events.  The  shape  of  the  particle  flux  decline  is  of  particular 
importance:  power-law  time  dependence  indicates  the  dominance  of  diffusive 
propagation,  whereas  exponential-law  decline  emphasizes  convection  transport. 
Depending  on  the  solar  wind  speed  both  exponential  (at  <10  MeV)  and  power-law  (>30 
MeV)  declines  can  be  present  in  the  same  event.  A  statistical  investigation  of  SEP 
events  extended  for  a  long  period  suggests  that  about  90  %  of  SEP  decays  are 
characterised  by  exponential  declines.  Distributions  of  the  total  durations,  T,  of 
exponential  declines  with  and  without  shocks  differ  insignificantly  for  T>10  hrs.  Values 
of  the  decay  time  obtained  theoretically  are  reasonably  close  to  the  fitted  slopes  in  nearly 
half  of  all  cases  if  one  uses  the  average  V  solar  wind  speed  values  measured  when  the 
corresponding  plasma  in  which  particles  were  convected  arrives  to  the  observer. 
Dependences  as  the  variation  of  the  decay  time  with  energy  and  angular  distance  of  the 
observer  from  the  flare  heliolongitude  are  considered  as  well. 
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15.  Peroids  of  Quasi-Stationan’  Conditions  in  Interplanetarj  Space  According  to 
Sequences  of  SEP  Events 

E.l.  Daibog.  Yu. I.  Logachev.  S.  Kahler,  K.  Kecskemety 
Advanses  in  Space  Research,  2003,  in  press 

The  values  of  the  characteristic  decay  time  of  particle  fluxes  in  SEP  events  vary,  as  a 
rule,  considerably  from  event  to  event.  We  discovered,  however,  that  at  times  sequences 
of  events  having  similar  decay  times  were  observed  over  long  time  intervals  (up  to  one 
month,  and  even  longer  in  a  few  cases).  The  values  of  the  decay  times,  however,  differed 
in  different  sequences.  Such  sets  of  events  took  place  during  nearly  all  superevents  in 
the  2 1  st  cycle  of  solar  activity.  However,  we  can  see  a  few  time  intervals  characterized 
by  similar  flux  declines  in  periods  between  the  superevents  as  well.  Thus,  it  turns  out 
periodically  that  either  the  interplanetary  medium  parameters  determining  conditions  of 
energetic  particle  propagation  (diffusion,  convection  and  adiabatic  deceleration)  or  the 
conditions  of  trapping  of  particles,  accelerated  in  flares  and  at  interplanetary  shocks, 
between  strong  magnetic  fields  near  the  Sun  and  shock  front,  stay  quite  stable  for  long 
periods  of  time.  The  investigation  of  such  sequences  of  events  can  also  be  useful  for 
forecasting  characteristics  of  SEP  events. 


16.  Statistical  Study  of  Spectral  Characteristics  of  SEP  Event  Flux  Declines  during 
1974-2001. 

E.l.  Daibog,  Yu.I.  Logachev,  S.  Kahler,  K.  Kecskemety,  S.  McKenna-Lawlor 
Advanses  in  Space  Research,  2003,  in  press 

The  interplanetary  space  is  not  a  passive  medium  for  propagating  previously  accelerated 
particles,  but  in  itself  influences  the  particles  by  changing  their  energies  due  to 
interactions  with  magnetic  field  inhomogeneities.  These  processes  are  manifested  in  the 
energy  spectra  of  SEP  events.  Here  the  integral  energy  fluxes  of  protons  in  the  energy 
range  4-30  MeV  are  investigated  on  the  basis  of  two  data  sets.  The  first  is  a 
homogeneous  set  of  data  obtained  by  the  CPME  instrument  aboard  the  IMP  8  satellite 
between  1974  and  2001  for  all  SEP  events  where  >10  MeV  proton  fluxes  exceed  10 
particle/cm^  s  sr.  The  second  consists  of  lower  energy  proton  fluxes  in  the  range  of  0.75- 
6  MeV  measured  by  the  LION  instrument  aboard  SOHO  between  1995  and  2001.  Other 
important  characteristics  of  SEP  events  include  the  rates  of  decrease  of  particle  flux, 
which,  as  well  as  peak  flux  time,  is  an  integral  feature  of  the  interplanetary  medium 
within  a  considerable  region,  surrounding  the  observation  point.  It  is  shovm  that  the 
average  values  of  energy  spectral  index,  y,  and  characteristic  decay  time,  x,  are  all 
changing  with  the  solar  activity  phase.  Distributions  of  y  and  x  are  obtained  in  SEPs  with 
and  without  shocks  (and  during  different  phases  of  events:  just  after  peak  flux  and  late 
after  maximum). 
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List  of  presentations  at  conferences  and  meetings  with  abstracts 


1.  27-th  International  Cosmic  Ray  Conference,  Hamburg,  August  2001 

Energetic  Electron  Spectra  in  Solar  Energetic  Particle  Events  Associated  and 
Non-Associated  to  Coronal  Mass  Ejections 

V,  Stolpovskii,  E.  Daibog,  G.  Erdos,  S.  Kahler ,  K.  Kecskemety  and  H.  Kunow 

We  considered  energetic  electron  spectra  in  solar  energetic  particle  (SEP)  events  after 
flare/coronal  mass  ejection  (CME)  associations  and  after  flares  only  using  measurements 
onboard  various  spacecraft  (s/c)  at  radial  distances  of  0.3-1  AU  during  the  21st  solar 
activity  cycle.  Statistics  included  about  a  hundred  events  of  both  types.  More  than  50 
events  of  each  type  were  observed  simultaneously  at  various  points  of  the  inner 
heliosphere.  Energy  spectra  in  the  range  0.3  to  3  MeV  were  generated  from  maximum 
flux  at  each  energy.  The  relationship  between  the  exponent  of  integral  electron 
spectrum,  y  and  CME  speed,  V  and  angular  width  was  considered  and  it  was  obtained 
that  electron  spectra  become  harder  with  increase  of  V.  The  best  correlation  between  y 
and  V  arrives  when  the  observer's  magnetic  footpoints  are  in  the  limits  of  the  CME 
angular  width.  In  this  case  the  best  fit  approximation  of  y(V)  looks  as  yocV°’^  with  y  and 
V  in  the  range  of  4.5-1 .5  and  500-2000  km/s,  correspondingly. 


Decay  Phases  in  Gradual  and  Impusive  Solar  Energetic  Particle  Events 

E.  Daibog,  V.  Stolpovskii,  G.  Erdos,  S.  Kahler,  K. Kecskemety  and  H.  Kunow 

We  have  studied  the  decay  phases  of  energetic  particle  (SEP)  events  associated  and  non- 
associated  to  coronal  mass  ejections  (CMEs)  and  interplanetary  (IP)  shocks  using 
multispacecraft  observations  in  the  end  of  1970s  and  1980s.  Here  several  examples  of 
successive  SEP  events  observed  simultaneously  at  widely  spaced  s/c  are  presented.  It  is 
pointed  out  that  often  during  long  time 

intervals  (up  to  one  month)  the  decays  of  >4  MeV  proton  and  >0.3  MeV  electron 
intensities  are  very  similar  over  a  wide  range  of  angles  between  the  observer  and 
proposed  source  (up  to  60  o  and  more).  The  spatial  and  temporal  invariance  of  the 
particle  intensities  during  the  late  phase  of  large  proton  events  associated  to  the  fastest 
CMEs  and  CME-driven  shocks  is  valid  not  only  in  a  wide  angular 
regions  but  at  different  distances  from  the  Sun  from  0.3  to  1  AU. 
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2.  18-th  Eurorcan  Cosmic  Ray  Symposium,  Moscow,  July  2002 


Statistical  characteristics  of  SEP  event  flux  declines  through  a  long  period  of  time 
(1974-2001) 

Daibog  E.I.,  Logachev  Yu. I.,  Kahler  S.W.,  Kecskemety  K. 

The  rate  of  decrease  of  particle  flux,  as  well  as  peak  flux  time,  is  an  integral  character  of 
the  interplanetary  medium  within  a  considerable  volume,  surrounding  an  observation 
point.  We  performed  a  statistical  analysis  of  the  rates  of  decrease  in  all  SEP  events  with 
>4  MeV  proton  flux,  exceeding  2-3  part/cm^  s  sr,  on  the  basis  of  a  homogeneous  data  set 
obtained  by  the  CPME  instrument  aboard  the  IMP  8  satellite  between  1974  and  2001. 
We  consider  the  intervals  of  intensity  time  profiles,  during  which  declines  were 
characterized  by  exponential  law  typical  for  <50  MeV  protons  of  most  SEP  events.  The 
distribution  of  characteristic  decay  times,  x,  of  >4  MeV  and  >10  MeV  protons  for  355 
events  has  been  built. 

The  comparison  of  a  complete  set  of  events  from  21  and  22  solar  activity  cycles 
indicates  that  both  the  total  number  of  events  and  the  average  values  of  x  differed  in 
these  cycles.  Integral  energy  spectra  of  proton  fluxes  in  the  energy  range  4-30  MeV  were 
also  obtained  for  the  above  time  periods.  In  addition,  we  analyze  the  correlation  of  x,  the 
durations  of  exponential  declines  and  spectral  exponents  with  interplanetary  medium 
conditions,  in  particular  with  CMEs  and  interplanetary  shocks. 


3.  COSPAR,  Houston,  USA,  August  2002 


Peroids  of  Quasi-Stationary  Conditions  in  Interplanetary  Space  According  to 
Sequences  of  SEP  Events 

E.I.  Daibog,  Yu.I.  Logachev,  S.  Kahler,  K.  Kecskemety 

The  values  of  the  characteristic  decay  time  of  particle  fluxes  in  SEP  events  vary,  as  a 
rule,  considerably  from  event  to  event.  We  discovered,  however,  that  at  times  sequences 
of  events  having  similar  decay  times  were  observed  over  long  time  intervals  (up  to  one 
month,  and  even  longer  in  a  few  cases).  The  values  of  the  decay  times,  however,  differed 
in  different  sequences.  Such  sets  of  events  took  place  during  nearly  all  superevents  in 
the  21st  cycle  of  solar  activity.  However,  we  can  see  a  few  time  intervals  characterized 
by  similar  flux  declines  in  periods  between  the  superevents  as  well.  Thus,  it  turns  out 
periodically  that  either  the  interplanetary  medium  parameters  determining  conditions  of 
energetic  particle  propagation  (diffusion,  convection  and  adiabatic  deceleration)  or  the 
conditions  of  trapping  of  particles,  accelerated  in  flares  and  at  interplanetary  shocks, 
between  strong  magnetic  fields  near  the  Sun  and  shock  front,  stay  quite  stable  for  long 
periods  of  time.  The  investigation  of  such  sequences  of  events  can  also  be  useful  for 
forecasting  characteristics  of  SEP  events. 
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Statistical  Study  of  Spectral  Characteristics  of  SEP  Event  Flux  Declines  during 
1974-2001. 

E.I.  Daibog.  Yu. I.  Logachev,  S.  Kahler,  K.  Kec.skemety,  S.  McKenna-Lawlor 

The  interplanetary  space  is  not  a  passive  medium  for  propagating  previously  accelerated 
particles,  but  in  itself  influences  the  particles  by  changing  their  energies  due  to 
interactions  with  magnetic  field  inhomogeneities.  These  processes  are  manifested  in  the 
energy  spectra  of  SEP  events.  Here  the  integral  energy  fluxes  of  protons  in  the  energy- 
range  4-30  MeV  are  investigated  on  the  basis  of  two  data  sets.  The  first  is  a 
homogeneous  set  of  data  obtained  by  the  CPME  instrument  aboard  the  IMP  8  satellite 
between  1974  and  2001  for  all  SEP  events  where  >10  MeV  proton  fluxes  exceed  10 
particle/cm  s  sr.  The  second  consists  of  lower  energy  proton  fluxes  in  the  range  of  0.75- 
6  MeV  measured  by  the  LION  instrument  aboard  SOHO  between  1995  and  2001.  Other 
important  characteristics  of  SEP  events  include  the  rates  of  decrease  of  particle  flux, 
which,  as  well  as  peak  flux  time,  is  an  integral  feature  of  the  interplanetary  medium 
within  a  considerable  region,  surrounding  the  observation  point.  It  is  shown  that  the 
average  values  of  energy  speetral  index,  y,  and  charaeteristic  decay  time,  t,  are  all 
changing  with  the  solar  activity  phase.  Distributions  of  y  and  x  are  obtained  in  SEPs  with 
and  without  shocks  (and  during  different  phases  of  events:  just  after  peak  flux  and  late 
after  maximum). 


4.  Pisa  Solar  Wind  10  conference,  Pisa,  Italy,  October  2002 

Invariance  of  proton  time  profiles  during  the  late  phase  of  SEP  events: 
multispasecraft  observations. 

Daibog  E.I.,  Stolpovskii  V.G.,  Kecskemety  K. 

Using  simultaneous  observations  from  Helios-1, -2  and  IMP-8  s/c  we  examined  the 
properties  of  spatial  and  temporal  invariance  of  Ep>4  MeV  proton  and  Ee>0.3  MeV 
electron  spectra  during  the  decay  phase  of  SEP  events  after  flare-CME  -  CME-drived 
shock  associations.  The  spectra  and  decay  times  of  charged  particles  are  nearly  invariant 
over  large  angular  distances  from  the  flare  point.  According  to  the  west  or  east  position 
of  s/c  relative  to  shock  nose  we  observe  delay  or  advance  of  invariance  onset  times 
following  the  shock  arrival  times  increasing  with  angular  distance  growth. 
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5.  EGS-AGU-EUG  Joint  Assembly  Geophysical  Research),  Nice,  France,  April 
2003 

Diagnostic  of  the  state  of  the  interplanetary  medium  according  to  the  decay  phase 
of  SEP  events. 

Daibog  E.I.,  Kecskemety  K. ,  Kahler  S.,  Logachev  Yu.I.  EGS-AGU-EUG  Joint 

Generalized  parameters  characterizing  the  interplanetary  medium  (IM)  state  are  the 
shape  and  the  rate  of  decline  of  charged  particle  flux  in  solar  energetic  particle  (SEP) 
events  (just  as  the  temperature  of  persons'  body  to  some  extent  is  a  generalized 
characteristics  of  his  state  of  health). 

The  shape  of  the  particle  flux  decline  is  of  particular  importance:  power-law  time 
dependence  says  about  prevalence  of  diffusive  propagation  of  particles,  exponential-law 
decline  about  convective  processes.  Depending  on  the  solar  wind  speed  both  kinds  of 
declines  may  be  realized  in  the  same  event:  exponential  for  low  energy  (<  10  MeV) 
particles  and  power-law  for  high  energy  (>30  MeV)  ones. 

Statistical  investigation  of  SEP  events  during  a  long  period  of  time  (from  1968  to  2001) 
shows  that  90  %  of  SEP  decays  are  characterised  by  exponential  declines.Distributions 
of  the  values  of  characteristic  decay  time,  x,  shows  that  most  part  of  the  time  during 
which  diagnostics  is  possible  the  interplanetary  medium  is  in  a  "basically  disturbed" 
state  characterised  by  t=16  ±  6  hrs.  This  permits  to  say  about  predominant  quasi  stable 
state  of  the  interplanetary  medium,  which  takes  about  one  half  of  the  total  time. 

As  SEP  events  usually  last  few  days,  the  rest  time  such  a  diagnostics  is  inaccessible. 
However,  often  SEPs  follow  one  by  one  and  all  events  in  the  sequences  have  similar 
decay  times,  x,  which  indicates  stability  of  conditions  in  the  IM  during  the  hole 
sequence.  This  permits  to  analize  the  IM  state  not  only  for  2-4  days  of  individual  event 
but  during  2-3  weeks  and  longer.  The  most  prolonged  sequence  and  IM  stable  state  was 
observed  in  1972  by  IMP-5  and  Prognoz-1.  Duration  of  this  period  took  above  2 
months,  i.e.  more  than  2  full  Solar  rotations. 


6.  International  Solar  Cycle  Studies  Symposium  2003:  "Solar  variability  as  an 
input  to  the  Earth's  environment,  June23-28,  Tatranska  Lomnica,  Slovak 
Republic 

Statistical  Properties  of  the  Decay  Phase  of  SEP-events 

Kecskemety  K.,  Daibog  E.I.,  Kahler  S.W.,  Logachev  Yu.I. 

Generalized  parameters  characterizing  the  slate  of  the  interplanetary  medium  (IM) 
include  the  functional  form  and  the  rate  of  decline  of  charged  particle  flux  in  solar 
energetic  particle  (SEP)  events.  The  shape  of  the  particle  flux  decline  is  of  particular 
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importance:  power-law  time  dependence  indicates  the  dominance  of  diffusive 
propagation,  whereas  exponential-law  decline  emphasizes  convection  transport. 
Depending  on  the  solar  wind  speed  both  exponential  (at  <10  MeV)  and  power-law  (>30 
MeV)  declines  can  be  present  in  the  same  event.  A  statistical  investigation  of  SEP 
events  extended  for  a  long  period  suggests  that  about  90\%  of  SEP  decays  are 
characterised  by  exponential  declines.  Distributions  of  the  total  durations,  T,  of 
exponential  declines  with  and  without  shocks  differ  insignificantly  for  T>10  hrs.  Values 
of  the  decay  time  obtained  theoretically  are  reasonably  close  to  the  fitted  slopes  in  nearly 
half  of  all  cases  if  one  uses  the  average  V  solar  wind  speed  values  measured  when  the 
corresponding  plasma  in  which  particles  were  convected  arrives  to  the  observer. 
Dependences  as  the  variation  of  the  decay  time  with  energy  and  angular  distance  of  the 
observer  from  the  flare  heliolongitude  are  considered  as  well. 
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Attachments 


Table  lA  p.  64 

Table  2 A  p.  76 
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Table  lA 


Table  of  solar  cosmic  ray  events  according  to  the  data  of  CPME 
instrument  from  the  IMP-8  spacecraft  during  1974-2001 


The  columns  are: 

1.  The  date  of  the  beginning  of  the  time  interval  (typically  10-15  days),  covering  a  group 
of  events. 

2.  Beginning  of  the  exponential  period  during  the  event 

3.  Jmax  -  proton  Ep  >  10  MeV  (cm^ssr) peak  flux  during  the  event. 

4.  Characteristic  decay  time  for  J(Ep>4  MeV)  for  the  exponential  period.  Index  “b”  - 
flux  at  the  end  of  period  AT  is  close  to  the  background. 

5.  The  same  as  4.  for  J(Ep>10  MeV). 

6.  Exponent  of  energetic  spectra  for  exponential  period. 

7.  The  duration  of  the  exponential  period  . 

8.  The  marks  are 

P  -  the  first  exponential  period  after  maximum  of  the  event 
W  -  exponential  period  without  shocks 
Sh  -  exponential  period  with  the  shock 

I  -  exponential  period  with  the  shock  not  influenced  on  the  rate  of  intensity 
decline. 


The  values  oft ,  Ydecay  and  AT  were  obtained  from  the  time  profiles  of  proton  fluxes 

Table  lA  includes  information  on  208  SEP  events  during  the  period  from  1974  to  2001, 
with  the  peak  fluxes  exceeded  2-3  protons/cm^.s.ster. 


1 

2 

3 

4 

5 

6 

7 

8 

Years,  date 

Month, 

day 

Jmax 

(>10MeV) 

T>  4 
MeV 

T  >  10 
MeV 

Tdecay 

AT, 

hours 

Note 

1974  06  05 

06.  09 

3 

2,1 

3,6 

4,4 

1 

P,  W 

06.  10 

— 

2,0 

8,9 

2,1 

2 

Sh 

1974  06  26 

06.29 

0,2 

9,6 

29,2 

1,7 

18 

P,  W 

07.04 

30 

7,5 

8,1 

2,4 

12 

P,  Sh,  I 

07.  06 

300 

2,6 

4,4 

2,1 

2 

P,  Sh 

07.  06a 

— 

17,2 

14,3 

1,9 

36 

W 

07.08 

12,0 

15,1 

2,1 

40 

— 

1974  08  18 

08.21 

0,4 

25,0 

78,0 

1,1 

20 

W 

1974  09  10 

09.  14 

100 

27,0 

18,2 

2,0 

30 

P,  w 

09.  16 

2 

5,2 

7,8 

2,0 

12 

Sh 

09.  18 

1 

14,5 

24,8 

1,1 

24 

— 

09.21 

100 

15,9 

12,5 

2,3 

27 

P,  w 

64 


09.21a 


09.  25 


09.  25a 


09.  27 


5.7  2,1 


23.4  0,6 


.8 


49,7  0.6 


3,5  1,2 


1,1 


1,6 


1,6 


2,0 


8,8  I  0,7 


26,0 

22,0 

26,0 

52,0  b 

14,4 

37,8 

37,5 

81,0  b 

16,2 

21,6 

16,2 

21,6 

44,0  b 

150  b 

8,3 

11,1 

65 


66 


1980  01  10 


04.  02 


iraE 

■ 

19.2 

2.0 

20 

P.  W 

8J 

■E 

1,1 

20 

P,  W 

68 


mmmm 


06.06 


06.  09 


06.  12 


I  06.  15 


06.  16 


06.  20 


06.  22 


1982  06  23  06.30 


1982  07  09  07.13 


07.  14 


07.  23 


07.  23a 


07.  26 


09.  06 


09.  08 


8.0  I  27,0  2,1 


1982  07  21 


1982  09  03 


WfeWlIj 


1982  11  23 


1982  1205 


1982  12  13 


1982  12  23 


1983  01  01 


1983  02  01 


1983  03  11 


1983  04  14 


1983  05  09 


1983  06  14 


09.22 


11.24 


11.24a 


11.25 


11.27 


11.27a 


11.30 


12.  08 


12.  09 


12.  12 


12.  18 


12.  19 


12.20 


12.  22 


12.  28 


12.  30 


01.02 


01.06 


01.09 


02.  04 


02.  07 


03.  14 


04.  16 


04.  25 


05.  12 


05.  13 


05.  14 


05.  16 


05.  17 


06.  15 


06.  17 


06.  24 


06.  28 


115 

115 

1,2 

50 

22,6 

39,8  b 

1,8 

24 

25,0  b 


,5 


22,0 


3,9 


15,0 


8,3 


1,7 


4,7 


0,8 


24,3  b 


2,5 


2,1 


33,0 


43,0 


24,1 


8,1 


5,0 


22,6 


12,5 


12,2 


43,5 


11,7 


10,0 


47,0  b 


8,9 


13,0 


70,0 


34,0  b  0,8 


2,1 


60.0  b  1,1 


29,0 


36,0  b  I  35,0  b 


10,4 


22,0  10,4 


7,8 


15,1 


18,3  b 


8,3 


40,0  b 


31,3  b 


12,2 


55,0  b 


19,0  b 


22,4  b 


78,0  b 


40,0  b 


37,5 


3,1 


104  b 


24,0 

36 

IKSl 

1,9 

8 

18,7 

1,7 

24 

11,2 

1,7  1 

30 

2 

12.2 

3.3 

10 

69 


70 


1989  03  07 

03.08 

03.  10 

1989  03  17 


1989  04  10 


1989  05  01 


1989  05  20 


1989  06  15 


1989  06  25 


1989  07  21 


1989  08  12 


1989  09  28 


1989  10  19 


1989  11  03 


1989  11  15 


03.  13 


03.  14 


03.  19 


03.20 


03.24 


03.  25 


04.  12 


04.  15 


05.  05 


05.  06 


05.07 


05.  09 


05.23 


05.24 


06.  18 


06.  19 


06.23 


06.  23a 


06.29 


06.  30 


07.01 


07.01a 


07.25 


08.  14 


08.  17 


08.20 


09.  05 


09.  06 


09.  13 


09.  14 


09.  18 


10.  03 


10.  20 


10.  24 


10.26 


10.29 


11.09 


11.  15 


11.  17 


11.20 


11.22 


11.28 


11.29 


12.01 


12.  04 


02.  03 


02.  03a 


11,4 

19,7 

1,6 

18 

6,8 

7,3 

0,9 

2 

15,6  b 

23,5  b 

0,9 

10 

P,  Sh 


P,  Sh,  I 


P,  W 


71 


1990  02  11  102.14 


P.  W 


P.  W 


P,  Sh 


P.  W.  I 


P,  Sh 


1991  06  15 


1991  06  27 


1991  07  07 


1991  08  25 


1991  09  03 


1991  09  29 


1991  1028 


1991  1221 


1992  02  06 


1992  02  25 


1992  03  07 


1992  05  08 


1992  06  25 


1992  08  04 


1992  10  30 


1992  11  29 


1993  03  04 


1994  02  19 


1994  10  18 


1995  10  19 


06.  13 


06.  16 


06.  18 


07.  03 


07.  09 


07.  10 


07.  12 


08.28 


08.  28a 


09.  07 


09.  08 


10.  02a 


10.  28 


10.  29 


10.30 


10.  31 


12.  22 


12.  23 


12.  26 


12.29 


12.31 


01. 04 


01.  05 


02.  08 


02.26 


02.  26a 


02.27 


02.  28 


03.  08 


03.  11 


03.  16 


05.  10 


05.  11 


06.  27 


07. 


08.05 


08.07 


11.01 


11.05 


11.30 


03.05 


03.07 


03.  13 


02.21 


10.20 


10.21 


10.  20 


iim 

18,2 

1,3 

24,0 

17,0 

1,4 

13.5 

20.0 

1.45 

26,0 

lOHI 

1.5 

9,7 

7,6 

7,6 

BBil 

6,8 

20,2 

22,9 

EEm 

11,3 

12,0 

6,4 

12,0  b 

2,1 

15,4 

17,5 

1,9 

23,2 

40,0 

2,1 

10,7 

28,2 

28,3  0,6 


,7 


10,4b  0,9 


13,2  I  1,6 


2,7 


— 

P.  — 

20 

P,  W 

26 

Sh 

36 

P,  W 

2 

P.  W 

16 

P,  W 

12 

P,  Sh 

8 

P,  Sh 

40 

10 

P,  W 

8 

P,  w 

20 

P,  Sh 

28 

2 

P,  w 

22 

p,  -- 

P,  W 


P,  Sh 


P,  W 


P,  W 


P,  Sh 


P,  W 


P,  W 


P,  W 


P,  W 


P,  W 


1,25 

50 

1,6 

18 

1,4 

44 

73 


11,1 


24,0 


19,2 


29,5  b 


1,4 


1,7 


16,2 


1,8 


27,0 


12,0 


55,0  b  1,8 


7,3  b  1,3 


39,5  b  2,1 


2,25 


80 

P,  W 

18 

P,  W 

2 

W 

3 

P,  Sh 

30 

P,  W 

2 

P,  Sh 

26 

P,  W 

11.21 

.... 

40,0  b 

— — 

.... 

1999  12  03 

12.  10 

0,9 

kkB 

22,5  b 

1,6 

8 

P,  W 

2000  04  02 

>20 

3.9 

4.5 

3.2  1  -- 

P.  W 

5 

11,5 

12.5 

3.2 

6 

2000  06  04 

06.09 

40 

10,4 

12,2 

3,8 

14 

P,  w 

06.  11 

30 

12,5 

11,2 

1.6 

8 

P,  Sh 

06.  12 

— — 

8,1 

15,6  b 

1,6 

8 

P,  -- 

2000  07  12 

07.  13 

3 

29,1 

2,9 

12 

P,  Sh 

07.  16 

200 

12,5 

14,1 

1,6 

20 

P,  Sh 

07.  17 

— - 

27,8 

26,0 

1.6 

26 

— 

07.23 

10 

kh 

8,4  b 

1,6 

12 

P,  W 

2000  09  1 1 

09.  14 

150 

17,8 

18,2 

mm 

24 

09.  17 

10 

3,9 

6,0 

““““ 

09.  18 

92,0 

24 

Sh,I 

09.  20 

1 

24,8  b 

52,0  b 

1,2 

18 

P,  W 

2000  10  15 

10.  17 

10 

mm 

23,0 

1,6 

70 

P,  W 

10.  27 

8 

17,5 

27,2 

2.7 

26 

P,  W 

2000  1 1  07 

11.09 

10000 

7,3 

1.6 

27 

P,  W 

11.  10 

500 

— 

Sh 

11.  11 

— 

1,7 

2,1 

DH 

8 

— 

2000  1 1  23 

11.25 

80 

mm 

1,5 

8 

P,  W 

11.27 

800 

20,9 

16,2 

1,5 

36 

P,  Sh,  I 

2001  01  26 

01.30 

30 

8,3 

48 

01.31 

2 

8,9  b 

1,6 

Sh 

2001  03  28 

03.31 

20 

■ssai 

2,0 

P,  Sh 

04.  05 

100 

22,2 

50 

P,  Sh 

2001  04  09 

04.  12 

100 

5,0 

04.  13 

40 

13,8 

KEBI 

1,3 

44 

P,  Sh 

04.  16 

800 

wmM 

10,0 

1,5 

48 

P,  W 

04.  18 

200 

11,1 

10,0 

lOH 

24 

2001  05  18 

05.21 

>8 

20 

15 

— 

2001  06  10 

06.  16 

20 

16,0 

13,9 

2,1 

16 

P,  W 

On  the  basis  of  this  table  distributions  of  different  parameters  of  SEP  events  were  built, 
allowing  to  obtain  statistically  well  grounded  conclusions.  Several  (from  1  to  7)  time 
intervals  were  distinguished  for  each  event  during  which  proton  flux  decayed  according 
to  the  exponential  low.  Each  of  them  was  considered  independently.  It  may  be  noted  that 
sometimes  for  long  lasting  events,  when  the  decay  phase  consisted  of  several 
exponential  periods,  it  was  possible  to  choose  one  longer  time  interval  (sometimes  it  was 
the  whole  decay)  which  could  be  described  by  a  single  characterstic  decay  time, 
interrupted  either  by  some  short  (forbush-like)  decrease  or  by  some  additional  flux, 
temporary  disturbing  general  exponential  behavior.  These  events  are  of  special  interest 
and  will  be  analyzed  in  detail  after  completion  the  statistical  study. 
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Flares  and  shocks  associated  with  SEP  events 
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